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Since the mid-1990s a new generation of Earth-observing satellites has been able to detect tropospheric 
air pollution at increasingly high spatial and temporal resolution. Most primary emitted species can be 
measured by one or more of the instruments. This review article addresses the question of how well we 
can relate the satellite measurements to quantification of primary emissions and what advances are 
needed to improve the usability of the measurements by U.S. air quality managers. Built on a compre- 
hensive literature review and comprising input by both satellite experts and emission inventory spe- 
cialists, the review identifies several targets that seem promising: large point sources of NO* and SO 2 , 
species that are difficult to measure by other means (NH 3 and CH 4 , for example), area sources that cannot 
easily be quantified by traditional bottom-up methods (such as unconventional oil and gas extraction, 
shipping, biomass burning, and biogenic sources), and the temporal variation of emissions (seasonal, 
diurnal, episodic). Techniques that enhance the usefulness of current retrievals (data assimilation, 
oversampling, multi-species retrievals, improved vertical profiles, etc.) are discussed. Finally, we point 
out the value of having new geostationary satellites like GEO-CAPE and TEMPO over North America that 
could provide measurements at high spatial (few km) and temporal (hourly) resolution. 
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1. Introduction 

In the past two decades major advances have occurred in the 
detection of atmospheric pollution from space. The generation of 
satellite instruments launched since 1995 has proved to be capable 
of observing a wide range of chemical species at increasingly high 
spatial and temporal resolution. In addition, the transformation of 
raw satellite retrievals to user-friendly, archived products has 
progressed considerably, such that the application of satellite ob- 
servations to a wide range of atmospheric problems is no longer a 
daunting prospect. Previous review papers (Borrell et al., 2003; 
Martin, 2008; Palmer, 2008; Hoff and Christopher, 2009) have 
stressed the value of space-based observations for measuring 
tropospheric air pollution from space and have pointed out the 
potential to infer emission strengths. Satellite instruments measure 
the concentration of a chemical species in a vertical column of air; 
and, if transport and chemical conversion do not greatly interfere or 
can be compensated for, then the measurements can be used to 
estimate the emission rate of the species directly under that col- 
umn of air. 

In 2011 the U.S. National Aeronautics and Space Administration 
(NASA) Applied Sciences Program (ASP) created an Air Quality 
Applied Sciences Team (AQAST) to serve the needs of U.S. air quality 
managers through the use of Earth science satellite data, suborbital 
data, and models (NASA, 2013a). One component of the AQAST 
activity is an assessment of the applicability to U.S. air quality 
management of research that links satellite retrievals with emis- 
sions. We pose two questions: How can air quality managers make 
use of satellite retrievals to improve emission estimates? and What 
developments are needed to enhance the usefulness of those re- 
trievals? It is not difficult to prepare a list of possible applications: 
verification of point-source compliance, quantification of uncertain 
area sources (including biogenic and other natural sources), esti- 
mation of neighboring emissions (say, in Canada and Mexico) and 
long-distance emissions (say, in China) that may influence U.S. air 
quality, the coordinated use of multiple-species retrievals, the 
detection of wildfires and other episodic events, etc. But in which of 
these areas does valuable capability lie? 

In support of U.S. air quality management, EPA has developed 
the National Emissions Inventory (NEI), a compilation of emission 
estimates collected from state, local, and tribal (S/L/T) air agencies, 
supplemented with other resources (EPA, 2013a). The NEI is 
released on a triennial cycle, with the most recent version for 2008. 
Some satellite data products are already included in the NEI, mainly 
related to wildfires. However, there are many other areas in which 
satellite products could potentially augment the NEI in the future. 

The purpose of this article is to review the literature and 
generate from it a U.S. perspective on what satellite observations 
can tell us about emissions. Because of this focus, the reader should 
not expect to find reference herein to the use of satellite data in 
studies of ambient concentrations of ozone or secondary aerosols 
or the movement of polluted air masses. Nor do we include studies 
outside North America unless they contribute significantly to 
relevant methods development or the novel treatment of emission 
sources that have the potential to be transferred to the U.S. context. 

The article is organized in five main sections. First, we review 12 
major satellite instruments that have the capability to detect and 
quantify emissions (Section 2). Next, we describe the derivation of 
retrievals from these instruments for eight important chemical 
species (Section 3). Methodological approaches to enhance source 
quantification are described in Section 4, together with the primary 
application studies that have been published to date, organized by 
source type. In Section 5, we draw together the satellite capability 
and the characteristics of the NEI to highlight promising avenues 
for future applications, which are summarized in Section 6. 


2. Satellites: platforms, instruments, and their characteristics 

In this section we briefly review the history of relevant satellite 
missions, highlighting those payload instruments that have helped 
to advance the study of tropospheric air pollution. We focus on 
modern instruments that have been widely used in the quantifi- 
cation of emissions or, at minimum, in the identification of source 
regions or the assessment of regional source strengths. We do not 
include older instruments or modern instruments that are of 
limited value for emissions estimation. 

On the European side, the modern era of pollution-observing 
satellites began with the launch by the European Space Agency 
(ESA) of the second European Remote Sensing satellite ERS-2 on 
April 21, 1995, that carried the Clobal Ozone Monitoring Experi- 
ment (GOME) designed to measure O3, N0 2 , and related species 
(Burrows et al., 1999). GOME was fully functional until 2003, after 
which only limited data retrieval was possible until 2011. The 
successor to the ERS mission, ENVISAT, was launched on March 1, 
2002, and carried the Scanning Imaging Absorption spectroMeter 
for Atmospheric CHartogaphY (SCIAMACHY) (Bovensmann et al., 
1999), which has been extremely valuable for observing a wide 
range of trace gases in the troposphere. Sadly, the contributions of 
SCIAMACHY ended when contact with ENVISAT was lost on April 8, 
2012. ESA and the European Organisation for the Exploitation of 
Meteorological Satellites (EUMETSAT) jointly initiated the MetOp 
mission for operational meteorology, and launched the MetOp-A 
satellite on October 19, 2006. Two instruments on MetOp-A are 
of relevance to this work: GOME-2 (Callies et al., 2000) — an update 
of the original ERS-2 instrument — and the Infrared Atmospheric 
Sounding Interferometer (IASI) (Clerbaux et al., 2009), which is 
capable of measuring CO, NH 3 , CH 4 , and other species. 

On the U.S. side, the modern era began on December 18, 1999, 
when NASA launched the Earth Observing System (EOS) satellite 
Terra, which featured three relevant onboard sensors: the Moder- 
ate Resolution Imaging Spectroradiometer (MODIS) (Barnes et al., 
1998), the Multi-angle Imaging SpectroRadiometer (MISR) (Diner 
et al., 1998), and the Measurements Of Pollution In The Tropo- 
sphere (MOPITT) instrument (Drummond, 1992). The first two of 
these instruments measure particulate matter (PM) and its optical 
effects, while the third specializes in the measurement of CO. Two 
and a half years later, NASA launched Aqua on May 4, 2002. Aqua 
carried a second MODIS instrument as well as the Atmospheric 
Infrared Sounder (AIRS) (Aumann et al., 2003). On July 15, 2004, 
NASA launched the Aura satellite to improve our understanding of 
the changing chemistry of Earth’s atmosphere. Two Aura in- 
struments are particularly valuable in an air pollution context: 
the Ozone Monitoring Instrument (OMI) (Levelt et al., 2006) and 
the Tropospheric Emission Spectrometer (TES) (Beer et al., 2001). 
Aqua and Aura are both components of a constellation of satellites 
called the “A-train” that makes coordinated measurements of the 
same location as they pass over it, one at a time, within an eight- 
minute period. Unfortunately, two subsequent satellites that 
would have further aided the investigation of tropospheric air 
pollution failed to reach orbit: the Orbiting Carbon Observatory 
(OCO) in February 2009 and Glory in March 2011. On January 23, 
2009, the Japan Aerospace Exploration Agency (JAXA) and collab- 
orating institutions launched the Greenhouse gases Observing 
SATellite (GOSAT) with the goal of measuring concentrations of C0 2 
and CH4 and their variability over time and space (Kuze et al., 
2009). Table 1 provides an overview of the major satellite plat- 
forms, their key instruments, and their primary characteristics. It 
shows the major species that each instrument is capable of 
observing. Further details are provided in subsequent sections of 
this paper to the extent that they are relevant to emissions 
quantification. 
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With no new launches of Earth-observing satellite instruments 
planned for the next few years, concern has been expressed that an 
interruption in the continuity of measurements could occur if 
present instruments were to fail. While this is undoubtedly true, 
prospects have improved in recent years with the approval of two 
new missions in the U.S. The next generation of Earth science sat- 
ellite missions was planned as part of the National Research Council 
2007 Decadal Survey. To study tropospheric trace gases and aero- 
sols from geostationary orbit, the NASA Geostationary Coastal and 
Air Pollution Events (GEO-CAPE) mission was recommended for 
launch in the 2013—2016 time frame (Fishman et al., 2012), which 
would measure SO 2 , NO 2 , and other ozone/aerosol precursor spe- 
cies hourly over North America with the highest ground resolution 
of 4 km x 4 km. Due to budget constraints, the GEO-CAPE launch 
was subsequently delayed beyond 2020 , though refinements to the 
mission plan might bring this date forward somewhat. Then, in 
November 2012, a new mission called TEMPO (Tropospheric 
Emissions: Monitoring of Pollution) was selected under NASA’s 
Earth Venture Instrument program with the specific goal of 
measuring tropospheric air pollutants (O 3 , NO 2 , SO 2 , formaldehyde, 
glyoxal, and aerosols) over North America with high spatial 
and temporal resolution from geostationary orbit (NASA, 2013b). 
TEMPO is planned to be launched as a hosted payload on a com- 
mercial satellite before 2020. Depending on priorities and budgets, 
TEMPO could be a precursor or component of the larger GEO-CAPE 
suite of instruments. 

For Europe, ESA is planning the Sentinel 4 and Sentinel 5 mis- 
sions, carried on EUMETSAT meteorological satellites, to measure 
atmospheric composition. The Sentinel 5 Precursor mission is 
planned to launch in 2015 with a TROPOMI instrument, which will 
measure NO 2 and SO 2 with improved ground resolution 
( ~ 50 km 2 ). Over Europe, the diurnal variations will be observed by 
the geostationary Sentinel 4 mission after 2018. Clerbaux and 
Crevoisier (2013) have expressed optimism about the future of 
lASI-type instruments for remote sensing in this timeframe. In 
addition, the Korea Aerospace Research Institute (KARI) is planning 
to launch an atmospheric chemistry satellite into geostationaiy 
orbit over East Asia in the 2018 time frame. 


3. Species observations 

This section describes the broad principles of satellite detection 
of atmospheric constituents, focusing on eight important emitted 
species that are observable from space at the present time. We 
review the conditions for measuring each species and some of the 
key studies undertaken to develop and apply the measurement 
methods. For the methods that are the furthest developed and the 
most used, we provide a little more detail on retrieval methods and 
archived products than for those that are primarily in research 
mode at the present time. The satellite instruments measure radi- 
ation over a wide region of the electromagnetic spectrum from the 
ultraviolet (UV) to the visible (Vis) in the range of ~ 0.23— 0.9 pm 
and in the infrared (IR) from ~ 4 to 50 pm. Instruments can be 
classified into two types, those that observe solar backscatter ra- 
diation in the UV— Vis and those that observe thermal IR emission. 
Chemical species in the atmosphere are detected by the absorption, 
or attenuation, of radiation of specific wavelengths along the path 
that the radiation travels through the atmosphere. The inferred 
concentration of the chemical species, known as a retrieval, is 
determined by a complex set of spectral fitting and radiative 
transfer calculations, well described by Martin (2008). For the N0 2 
retrievals described in Section 3.1, which are the most fully devel- 
oped and most widely used, we illustrate the process of deriving 
concentrations from radiation measurements in more detail, 
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leading to the production of operational retrievals, i.e., fully pro- 
cessed data sets available for user applications. 

3.1. Nitrogen dioxide 

Nitrogen oxides (NO*) are an important class of atmospheric 
species that are implicated in a number of environmental problems, 
including the formation of tropospheric ozone and aerosols, acid- 
ification, eutrophication, and human health effects. NO* is primarily 
formed by combustion processes, and emission factors vary ac- 
cording to the fuel type and combustion conditions. For this reason, 
NO* emissions can be quite uncertain, and there is potential for 
aiding the estimation of NO* emissions using space-borne in- 
struments. NO2 is the specific chemical form detected spectro- 
scopically, so the relationship between NO* emissions and N0 2 
concentrations must be taken into account (Lamsal et al., 2008, 
2011). Daytime global N0 2 datasets are produced daily using UV— 
Vis spectrometers aboard polar-orbiting, sun-synchronous satel- 
lites. OMI features the highest spatial resolution (13 km [along 
track] x 24 km [across track] at nadir) among this type of instru- 
ment and has been providing continuous global coverage since 
October 2004. The GOME-2 instrument has provided morning N0 2 
data since 2007. It is characterized by a spatial resolution of 
40 km x 80 km with global coverage in 1.5 days. In addition, the 
SCIAMACHY instrument produced N0 2 data from its launch in 2002 
until its demise in 2012. It was characterized by a spatial resolution 
of 30 km x 60 km but had coverage gaps due to alternating limb 
and nadir modes of observation, so global coverage was achieved 
every six days at the equator. Their predecessor instrument, GOME, 
provided N0 2 data from 1995 to 2003 at 40 km x 320 km resolu- 
tion, achieving global coverage every three days at the equator. 

Processing of satellite measurements typically involves three or 
more steps or processing “levels”. Calibrated Level 1 satellite data 
include high-resolution solar and earthshine spectra (typically 
better than 0.6 nm spectral resolution) for each individual mea- 
surement (pixel). Earthshine radiances are divided by a reference 
solar irradiance spectrum to linearize the N0 2 retrieval problem. 
The reflectance spectra are fitted to laboratory-measured trace-gas 
spectra (e.g., Vandaele et al., 1998) convolved with an instrumental 
slit function, a reference Ring spectrum (e.g., Chance and Spurr, 
1997), and a polynomial function that models the spectrally 
slowly varying scattering by clouds and aerosols and reflection by 
Earth’s surface. The fitting algorithms use the Differential Optical 
Absorption Spectroscopy (DOAS) method, in which a continuous 
light source is used to send and receive light through the atmo- 
sphere and measure narrow-band UV— Vis absorption features 
associated with trace-gas species. These are applied to a pre- 
selected “fitting window” in the spectral range of 0.4— 0.5 pm. 
The spectral fit yields a slant column density (SCD), which repre- 
sents the measured total areal density of N0 2 molecules along the 
effective optical path from the Sun through the atmosphere to the 
satellite. Additional processes include an air mass factor (AMF) 
calculation, which relates the trace gas abundance over the entire 
path length to the vertical column abundance, and partitioning of 
tropospheric and stratospheric components. Vertically resolved 
scattering weights (SW) or box air mass factors (box-AMF) describe 
the relationship between the slant column and the vertical column 
and are computed with a radiative transfer code that uses surface 
reflectivity, cloud fraction, cloud height, and surface topography. 
The profile-weighted average ofSW provides the total stratospheric 
or tropospheric AMFs used to convert measured SCD to the inferred 
Level 2 vertical column density (VCD). 

The overall error in the retrieval of tropospheric N0 2 columns 
arises from errors in the SCD, the stratosphere-troposphere sepa- 
ration, and the AMF calculation. The uncertainty in SCD and the 


stratosphere-troposphere separation dominates the total retrieval 
error over the oceans and remote areas. Over continental source 
regions, AMF errors arising from incorrect assumptions about the 
N0 2 vertical profile shape, surface reflectivity, cloud parameters, 
and aerosols dominate the overall retrieval errors. Estimated errors 
in the tropospheric N0 2 columns are ~30% under clear-sky con- 
ditions, but can reach up to 100% under certain conditions 
(Boersma et al., 2011; Bucsela et al„ 2013). Instrument-specific ar- 
tifacts, such as stripes in the OMI observations, may introduce 
additional errors (Celarier et al., 2008). The tropospheric N0 2 re- 
trievals and inferred emissions can benefit from the use of high 
spatial resolution information on surface reflectivity and N0 2 ver- 
tical profile shape in the AMF calculation, as used in the Berkeley 
High-Resolution (BEHR) retrieval (Russell et al., 2011). The influ- 
ence of spatial resolution on N0 2 retrievals has been investigated 
by Valin et al. (2011) and Heckel et al. (2011). 

Production of Level 3 N0 2 archived data typically includes re- 
gridding pixel data to a regular latitude-longitude grid, filtering 
the data, and time averaging to produce weekly-, monthly-, and 
seasonally-averaged data. Time averaging of the best-quality 
filtered data performed on a fine spatial grid (say, 1—5 km) allows 
better characterization of stationary point sources such as power 
plants or smelters. Pixel-averaging techniques result in “super- 
resolution” N0 2 maps, greatly enhancing the satellite detection 
capability for weak point sources (see Section 4.2). 

N0 2 column retrievals have been the subject of a number of 
validation studies involving intercomparison with a variety of 
ground-based measurements, aircraft measurements, and model 
outputs. These studies contain great detail about the conditions 
(seasons, locations, weather conditions, etc.) under which the sat- 
ellite retrievals can be expected to be the most reliable. These 
studies include Blond et al„ 2007; Boersma et al., 2008a; Bucsela 
et al., 2008; Hains et al., 2010; Huijnen et al., 2010; Lamsal et al., 
2008, 2010; and Martin et al., 2002, 2003. Different retrievals 
have been the subject of intercomparison (Boersma et al., 2008b; 
Lin et al., 2012). During the past decade, N0 2 retrieval algorithms 
have matured considerably, such that different algorithms now 
have good consistency (Boersma et al., 2011; Bucsela et al., 2013). 

3.2. Sulfur dioxide 

Sulfur dioxide (S0 2 ) is an atmospheric constituent that is linked 
to adverse respiratory effects and can also form secondary sulfate 
particles through chemical reactions in the atmosphere. Acid 
deposition is another detrimental effect of S0 2 . Satellite measure- 
ments that are sufficiently sensitive to anthropogenic S0 2 emis- 
sions are spectroscopically possible by measuring backscattered 
solar radiation in the UV. Measurements of anthropogenic S0 2 have 
been demonstrated using several of the instruments described in 
Section 2, including GOME (Eisinger and Burrows, 1998), GOME-2 
(Nowlan et al., 2011), SCIAMACHY (Lee et al., 2008b, 2009), and 
OMI (Krotkov et al., 2006, 2008; Yang et al., 2007). Currently, S0 2 
data from the OMI instrument are available in NASA Goddard public 
archives, such as GES-DISC and AVDC. 

Different algorithms have been developed to retrieve S0 2 from 
satellite data. DOAS S0 2 algorithms have to use spectral fitting 
windows at longer UV wavelengths (>315 nm), where O3 and S0 2 
absorptions are weak (Lee et al„ 2008b, 2009). Alternatively, the 
operational OMI S0 2 band residual difference (BRD) algorithm 
(Krotkov et al., 2006) uses shorter wavelengths between 310.8 nm 
and 315 nm that are most sensitive to S0 2 , explicitly accounting for 
the strong O3 absorption in this spectral region. In the absence of 
information concerning the vertical distribution of S0 2 , a constant 
AMF value of 0.36 is assumed in the Level 2 OMS02 PBL product. This 
value is based on an average summertime S0 2 profile measured by 


D.C. Streets et al. / Atmospheric Environment 77 (2013) 1 011-1042 


1015 


aircraft over the U.S. Mid-Atlantic region (Taubman et al., 2006). 
Although advanced SO 2 algorithms such as optimal estimation 
(Nowlan et al., 2011) and direct spectral fitting (Yang et al., 2010) 
have been demonstrated, they are not yet used operationally. 

Several factors influence the quality of the satellite-retrieved SO 2 
product. The sensitivity of satellite retrievals to SO 2 in the boundary 
layer is low, due to enhanced Rayleigh scattering, large O 3 absorp- 
tion, and low surface albedo (0.01—0.03, except for snow-covered 
surfaces). Level 1 radiance spectra (<320 nm) have lower signal- 
to-noise ratio and higher stray light compared to measurements at 
longer UV— Vis wavelengths (van den Oord et al., 2006). As a result, 
the noise of the SO 2 retrieval is high (one standard deviation = 1 — 
2 DU) even under optimal viewing conditions (i.e., no clouds, near 
nadir, low solar zenith angle, mid- to low-latitudes). Spatial aver- 
aging of pixel data reduces the noise, although not as effectively as 
would be expected for random noise (Krotkov et al., 2008). In the 
OMI Level 2 OMS02 data, there is no correction for clouds, so cloud 
screening and other filters (e.g., high solar zenith angle, snow cover, 
row anomaly, etc.) — as well as various adjustments such as correc- 
tions for the Pacific sector and local AMF (Lee et al., 2009, 2011a) and 
local bias (Fioletov et al., 2011) — are recommended for users. All 
necessary filters and corrections are currently applied in OMI Level 3 
processing, and users are encouraged to use the Level 3 daily gridded 
(0.25° x 0.25°) S0 2 product (OMS02e). 

Despite these limitations, the operational OMI S0 2 data have 
been utilized in a wide range of studies (Krueger et al., 2009), in 
combination with aircraft MAX-DOAS measurements (Walter et al., 
2012), in situ measurements (Dickerson et al., 2007; Lee et al., 
2011b), and chemical transport modeling (CTM) (Lee et al., 2011a). 
Topics studied include the transport and transformation of emitted 
SO 2 (Dickerson et al., 2007; Lee et al., 2008a; Li et al., 2010a; de Foy 
et al., 2009), top-down constraints on emissions and lifetime (Lee 
et al., 2008a, 2011a); spatiotemporal variations (McLinden et al., 
2012; Jiang et al., 2012), and the impacts of anthropogenic emis- 
sion controls on air pollution (Witte et al., 2009; Li et al., 2010b; 
Fioletov et al., 2011). Temporal averaging over a month or longer 
can effectively reduce the uncertainties in the daily S0 2 operational 
data, which are estimated at ~100% (Krotkov et al., 2008) and 
considered to be too large for most air quality applications. Using 
pixel averaging and spatial filtration techniques, moderate point 
sources of S0 2 can be quantified (see Section 4.2). 

Spectroscopic detection of S0 2 in the infrared is also possible, 
though sensitivity is hindered by the strong absorption by water 
vapor in the same spectral region (Clarisse et al., 2008). Detection of 
ground-level S0 2 by IASI has been demonstrated in principle, but 
presently the sensitivity is greatest at higher altitudes and for 
strong sources. Retrievals have been described by Clarisse et al. 
(2008) in the context of observations of the S0 2 plume from the 
2007 eruption of the Jebel at Tail' volcano in the Red Sea. 

3.3. Carbon monoxide 

CO retrievals have unique value for observing anthropogenic 
activity, evaluating emissions, and determining the transport and 
influence of these emissions downwind of sources. For these rea- 
sons, it has been one of the few mission-critical measurements in 
all aircraft campaigns of the NASA Global Tropospheric Chemistry 
Program (Fisher et al., 2010). CO is a U.S. EPA criteria pollutant, the 
sources of which fall into two main classes: secondary chemical 
production and incomplete combustion. Chemical production re- 
sults from hydrocarbon oxidation, while combustion sources 
include urban/industrial fossil-fuel burning, biofuel use, wildfires, 
and tropical biomass burning. CO is the dominant sink for the hy- 
droxyl radical (OH), and as such plays a critical role in controlling 
OH levels with implications for a wide range of atmospheric gases. 


CO is also a major precursor for tropospheric O 3 on a global scale. 
The CO global average lifetime against OH oxidation is about two 
months, providing sufficient time for CO to be an effective tracer of 
long-range pollution transport, yet sufficiently short for it to show 
large contrasts between polluted air masses and the background 
atmosphere (Edwards et al., 2004). 

Two CO spectral bands, the thermal infrared (TIR) fundamental 
at 4.6 pm and the near infrared (NIR) overtone band at 2.3 pm, can 
be used for passive remote sensing (Edwards et al., 1999). In- 
struments on each of NASA’s EOS platforms make nadir measure- 
ments of CO at 4.6 pm from low Earth orbit: MOPITT (Deeter et al., 
2003), AIRS (McMillan et al., 2005), and TES (Rinsland et al„ 2006b). 
CO observations are also made by IASI on the MetOp-A satellite 
(George et al., 2009). These measurements are generally most 
sensitive to CO in the middle troposphere, with limited vertical 
profile information available only under favorable conditions (e.g., 
Deeter et al., 2004). Using solar reflection at 2.3 pm allows for 
retrieval of the total CO column with measurement sensitivity to all 
altitudes, including the lowermost troposphere. These have been 
made by SCIAMACHY (Buchwitz et al., 2006) and by MOPITT 
(Deeter et al., 2009). 

Currently, MOPITT is the only satellite instrument to demon- 
strate the instantaneous multispectral retrievals necessary for an 
independent measurement of a trace-gas concentration in the 
lowermost troposphere (Worden et al., 2010). Retrieving CO from 
combined TIR and NIR radiances allows increased discrimination of 
vertical structure and upwards of three pieces of independent in- 
formation, as characterized by the degrees of freedom of signal 
(DFS) (Rodgers, 2000). The fact that significant information about 
CO concentration in the lowest 2 km is often available should 
improve emission estimates of not only CO, but also other com- 
bustion products for which emission ratios with CO can be devel- 
oped. Accurate emission inventories for CO are important for air 
quality management; yet, past atmospheric observations have 
often yielded results that were significantly at odds with emission 
inventories (Parrish, 2006; Hudman et al., 2008; Miller et al., 2008; 
Kopacz et al., 2010). Further results from CO satellite observations 
are presented in Section 4.5 on anthropogenic area sources, the 
source category where the strongest potential applications exist. 

3.4. Methane 

Methane (CH 4 ) is the second most important anthropogenic 
greenhouse gas after C0 2 , and, as a short-lived climate forcing 
agent with a lifetime of ~ 10 years, it provides a lever for slowing 
near-term climate change. Emission controls also have co-benefits 
for improving air quality (Fiore et al., 2002). Efforts to reduce CH 4 
emissions around the globe are undertaken against the backdrop of 
poor quantitative understanding of CH 4 sources (Heimann, 2011). 
Major anthropogenic sources include natural gas extraction and 
use, coal mining, landfills, livestock, rice cultivation, and biomass 
burning. Wetlands are the largest natural source. Bloom et al. 
(2010) were able to estimate worldwide CH 4 emissions from wet- 
lands and rice paddies using SCIAMACHY coupled with other 
measurements. Present-day global emissions of CH 4 are 
580 ± 100 Tg yr ', constrained by knowledge of the main global 
sink (reaction with the OH radical) from the methylchloroform 
budget (IPCC, 2007). However, the allocation of global emissions to 
source types and source regions is highly uncertain (Dlugokencky 
et al., 2011). To cite just one contemporary example, U.S. natural 
gas production has boomed in recent years as a result of new 
drilling techniques to unlock fuel trapped in shale formations. 
There is poor understanding of the associated CH 4 emissions, and 
the shifting nature of oil/gas exploration complicates their ac- 
counting in standard inventories. 
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Satellite observations provide a unique resource for constraining 
CH 4 emissions through inverse modeling, because of the dense and 
continuous data they provide. Methane can be retrieved from 
space-based measurements of solar backscatter in the shortwave 
infrared (SWIR) and of terrestrial radiation in the thermal infrared 
(T1R). SWIR retrievals are available from SCIAMACHY for 2003— 
2012 (Frankenberg et al., 2005, 2011), though pixel degradation in 
the detector reduced data quality after 2005, and from GOSAT for 
2009— present (Yokota et al., 2009; Butz et al., 2011; Parker et al., 
2011). The retrievals are sensitive to the entire tropospheric col- 
umn, but their dependence on reflected sunlight precludes obser- 
vations at night, over most ocean surfaces, and over most cloudy 
targets. TIR retrievals are available from AIRS for 2002— present 


(Xiong et al., 2009, 2010), TES for 2004— present (Worden et al., 
2012b), and IASI for 2007— present (Crevoisier et al., 2009; Razavi 
et al„ 2009); they have limited sensitivity to the lower tropo- 
sphere but can be performed at night, over land and ocean, and for 
cloudy scenes. 

SCIAMACHY data have been used in global and continental in- 
verse modeling studies of CH 4 emissions (Bergamaschi et al., 2007, 
2009; Meirink et al„ 2006, 2008), although the early studies were 
affected by retrieval biases propagating through the inversion 
(Bergamaschi et al., 2009). Fig. 1 compares SCIAMACHY observa- 
tions for June— August 2004 to concurrent boundary layer obser- 
vations from the NASA INTEX-A aircraft campaign, as well as to 
GEOS-Chem model results using a state-of-science emission 
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Fig. 1 . Methane concentrations over North America in June-August 2004 as observed by SCIAMACHY (column mixing ratios) and INTEX-A (mixing ratios below 700 hPa). The 
middle panels show the corresponding GEOS-Chem model values, and the bottom panels show the observed-model differences. 
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inventory as the a priori (Pickett-Heaps et al., 2011). There is 
excellent correlation between SCIAMACHY and INTEX-A vertical 
profile data (r = 0.93, regression slope = 0.96). Comparison with 
GEOS-Chem shows patterns of model underestimation in the cen- 
tral U.S. and overestimation along the east coast. Placing these 
discrepancies in the context of the a priori emissions suggests that 
livestock emissions are underestimated, while landfill and wetland 
emissions may be overestimated. But we also see that major 
sources in the south-central U.S. (likely from oil/gas recovery) 
appear to be missing from the emission inventory. Formal source 
attribution of the model overestimates and underestimates must 
be achieved through inverse modeling. Adjoint inverse modeling, 
as demonstrated by Meirink et al. (2006, 2008) and Bergamaschi 
et al. (2009), allows for emission estimates at the native resolu- 
tion of the model. Specific source-type attribution needs to be 
derived by comparing the high-resolution inverse emission esti- 
mates with detailed bottom-up information. 

The ability to validate satellite observations for CH 4 has greatly 
expanded in recent years with the increased availability of aircraft 
vertical profiles and of ground-based column observations from the 
Total Carbon Column Observing Network (TCCON). Wecht et al. 
(2012) used the extensive vertical profiles from the HIAPER Pole- 
to-Pole Observation (HIPPO) program over the Pacific to success- 
fully validate the new TES v5 retrieval, including information in the 
lower troposphere. Validation of GOSAT with TCCON data, 
including over North America, shows that it can capture column 
changes of a few tenths of a percent, providing strong constraints 
for inverse modeling (Butz et al., 2011 ; Parker et al., 2011 ). GOSAT is 
likely to be an important source for space-based observation of CH 4 
emissions in the future. 

The complementary vertical sensitivities of SWIR and T1R re- 
trievals can be combined to improve emission estimates by con- 
straining lower tropospheric CH 4 concentrations. Kuai et al. (2013) 
have demonstrated the utility of this approach using CO 2 retrievals 
from TES and TCCON, but to our knowledge no similar study has 
been done combining satellite-based CH 4 retrievals. Fig. 2 shows 
the flight paths of the CalNex aircraft campaign in April— June 2010 


and the locations of TES and GOSAT observations during this period. 
This figure is a good illustration of the relative spatial coverage of 
satellite and aircraft measurements that frame the intercompar- 
ison. By performing a set of inversions using combinations of Cal- 
Nex, TES, and GOSAT data during the CalNex period, it is possible to 
evaluate the utility of each instrument or combination of in- 
struments for constraining California CH 4 emissions. 

3.5. Nonmethane volatile organic compounds 

Nonmethane volatile organic compounds (NMVOC) in the at- 
mosphere have important implications for air quality and climate 
through their effects on atmospheric chemistry. They serve as 
precursors of organic aerosol, organic nitrates, and tropospheric O 3 , 
and can either enhance or deplete OH concentrations. Ability to 
observe NMVOC from space is mainly limited to Cj _3 molecules, 
because larger molecules do not have the fine band structure 
necessary for detection. Solar backscatter measurements of atmo- 
spheric columns require strong absorption features in the UV— Vis 
region and have been limited so far to formaldehyde (HCHO) 
(Chance et al., 2000) and glyoxal (CHOCHO) (Wittrock et al., 2006). 
Absorption bands in the 1R enable the measurement of many other 
species by solar occultation or terrestrial emission, including 
methanol (CH 3 OH), formic acid (HCOOH), hydrogen cyanide (HCN), 
ethane (C 2 H 6 ), ethylene (C 2 H 4 ), acetylene (C 2 H 2 ), acetonitrile 
(CH 3 CN), peroxyacetylnitrate (PAN, CH 3 C( 0 ) 00 NC> 2 ), and acetone 
(CH 3 C( 0 )CH 3 ). However, sensitivity for these measurements is 
limited to the middle or upper troposphere. 

Satellite observations of HCHO columns by solar backscatter are 
available from GOME, SCIAMACHY, OMI, and GOME-2. These in- 
struments provide a long and continuous record of observations 
(De Smedt et al., 2008; Stavrakou et al., 2009b). The detection limit 
is a few ppb in the boundary layer. HCHO is a high-yield product of 
NMVOC oxidation and has a lifetime of only a few hours against 
photolysis and oxidation. Satellite observations show large HCHO 
enhancements over vegetated and biomass burning regions, 
reflecting fast production from short-lived NMVOC. Urban areas are 



Fig. 2. Flight paths of the CalNex aircraft campaign in California in April-June 2010 and the locations of TES and GOSAT satellite observations during this period. 
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detectable only under very polluted conditions, reflecting the 
relatively small fluxes of short-lived anthropogenic NMVOC (Shim 
et al., 2005; Fu et al., 2007; Millet et al., 2008; Curd et al., 2010). 
Palmer et al. (2003) and Millet et al. (2006) have shown that HCHO 
in non-combustion regions is mainly a tracer of isoprene emissions 
because of its large flux, its short lifetime against oxidation by OH, 
and its prompt yield of HCHO as a first-generation product. 

A number of studies have demonstrated the value of HCHO 
columns from space to infer isoprene emissions in different parts of 
the world, using HCHO-isoprene relationships derived from a CTM 
(Palmer et al., 2003, 2006; Barkley et al„ 2008; Millet et al., 2008; 
Marais et al„ 2012). Although isoprene chemistry under low-NO x 
conditions is not well known, this appears to have little effect on 
HCHO yields (Palmer et al., 2006; Marais et al., 2012). The largest 
source of error appears to be the “smearing” associated with the 
time delay between isoprene emission and HCHO formation 
(Marais et al., 2012). This smearing error can be addressed with a 
formal inverse analysis using the adjoint of a CTM (Stavrakou et al., 
2009c). However, the results are then sensitive to the coupling 
between chemistry and transport on the time scales of boundary 
layer mixing. 

Glyoxal is measurable from space by solar backscatter in its 
420—460 nm absorption bands (Wittrock et al., 2006; Vrekoussis 
et al., 2009, 2010). Operational products are available from SCIA- 
MACHY and GOME-2. Like HCHO, it is produced in the atmosphere 
from the oxidation of NMVOC, and it has a lifetime of a few hours 
against photolysis and oxidation by OH. However, it has very 
different yields from HCHO, depending on the precursor. There is 
evidence for fast production of CHOCHO from the oxidation of 
isoprene, but the mechanism involved is not well understood. 
CHOCHO data from SCIAMACHY show prominent signatures from 
pyrogenic, biogenic, and anthropogenic sources, but also elevated 
concentrations over tropical oceans that suggest a marine biogenic 
source (Wittrock et al., 2006). Vrekoussis et al. (2009) examined the 
CHOCHO/HCHO column ratios observed by SCIAMACHY for 12 
continental source regions and found significant variability, 
reflecting different source signature mixes of precursor NMVOC. 

Several global model studies have focused on better defining the 
atmospheric budget of CHOCHO, including evaluation with in situ 
and SCIAMACHY observations (Fu et al., 2008; Myriokefalitakis 
et al., 2008; Stavrakou et al., 2009a). Models are unable to repro- 
duce the high values observed over tropical oceans and are also 
generally too low over land. Stavrakou et al. (2009a) argued that the 
latter underestimate could reflect a missing secondary biogenic 
source of CHOCHO accounting for 50% of the global CHOCHO 
source. Liu et al. (2012b) inferred a large missing source of CHOCHO 
over China from SCIAMACHY measurements, which they attributed 
to an underestimate of anthropogenic emissions of aromatics. 
Overall, the SCIAMACHY data imply shortcomings in our current 
understanding of CHOCHO sources, with potentially important 
implications for NMVOC emissions and organic aerosol formation. 

IR observations of NMVOC are available from the ACE-FTS 
(2003—) solar occultation instrument, the MIPAS (2002—) limb 
thermal emission instrument, and the TES (2004—) and IASI 
(2006—) nadir thermal emission instruments. Solar occultation and 
limb measurements are restricted to the upper troposphere. ACE- 
FTS data have been reported for CH 3 OH (Dufour et al., 2006), 
HCOOH (Rinsland et ai„ 2006a), C 2 H 4 (Herbin et al„ 2009), HCN 
(Rinsland et al., 2005), and C 2 H6 and C 2 H 2 (Rinsland et al., 2005; 
Gonzalez Abad et al., 2011). MIPAS data have been reported for 
C 2 H 6 (Clarmann et al., 2007; Glatthor et al., 2009), HCN (Glatthor 
et al„ 2009), HCOOH (Grutter et al., 2010), C 2 H 2 (Parker et al., 

2011) , CH 3 C(0)CH 3 (Moore et al„ 2012), and PAN (Wiegele et al., 

2012) . The ACE-FTS and MIPAS data provide climatological distri- 
butions and seasonal cycles in the upper troposphere, as well as 


correlations among species, that are of great value for testing global 
CTMs (Dufour et al., 2007; Gonzalez Abad et al., 2011). Their value 
for testing emission inventories is more limited at present because 
of the convolution with vertical transport and chemistry and 
because upper tropospheric observations integrate emissions in- 
formation over a large domain. Nevertheless, because of the diffi- 
culties in compiling bottom-up emission inventories of speciated 
NMVOC emissions from many source categories and in many parts 
of the world, the contribution of satellite observations of individual 
chemical species may ultimately prove to be of value. For example, 
Rinsland et al. (2007) used ACE-FTS data in freshly converted 
biomass burning plumes to infer emission ratios. 

Nadir IR observations from TES and IASI have higher horizontal 
resolution and are sensitive to lower altitudes, thus providing 
better constraints to relate observed concentrations to sources and 
to test emission inventories. TES has been used to detect large ur- 
ban enhancements of CH 3 OH (Beer et al., 2008) and to test current 
inventories of global biogenic CH 3 OH emissions (Wells et al., 2012). 
A new retrieval algorithm for TES observations of CH 3 OH presents 
seasonal and spatial variability (Cady-Pereira et al„ 2012). IASI has 
provided data for CH 3 OH, HCOOH, PAN, C 2 H 4 , C 3 H 6 , and furan 
(C 4 H 4 0) (Clarisse et al., 2011b). Detection of most of these species 
requires large enhancements from a point source such as a large fire 
(Clerbaux et al., 2009; Coheur et al„ 2009). However, high-quality 
global data sets have been produced for CH 3 OH and HCOOH 
(Razavi et al., 2011 ), and these have been used in inverse analyses to 
constrain emissions and secondary sources (Stavrakou et al., 
2011a, b). 

3.6. Ammonia 

Ammonia (NH 3 ) plays an important role in atmospheric chem- 
istry and biogeochemical cycles of reactive nitrogen (Galloway 
et al., 2003). By contributing to the transport and deposition of 
reactive nitrogen, NH 3 poses a hazard to sensitive ecosystems. The 
formation of ammonium sulfate and ammonium nitrate aerosol 
impacts human health, visibility, and climate. However, routine 
monitoring of surface-level NH 3 concentrations is sparse, as it is a 
difficult substance to measure in situ (von Bobrutzki et al., 2010). 
The potential for space-based observations of NH 3 is thus especially 
appealing. Profiles of atmospheric NH 3 have been retrieved from 
two instruments: TES (Beer et al., 2008; Shephard et al., 2011) and 
IASI (Clarisse et al., 2009, 2010), both of which measure thermal IR 
in the range of 940—969 cm -1 twice daily. 

The biggest challenges in remotely observing NH 3 are that much 
of the NH 3 is close to the surface, which makes it difficult to detect 
from space, and the horizontal distribution of NH 3 can have strong 
spatial gradients. Under conditions with significant thermal 
contrast near Earth’s surface and a large enough signal-to-noise 
ratio, both instruments can be sensitive to boundary-layer NH 3 
concentrations (Clarisse et al., 2010). The shape of the profiles 
retrieved by these instruments is strongly influenced by the 
retrieval a priori information, as the retrievals themselves typically 
contain less than one degree of freedom for signal. Typical retrieval 
sensitivities peak near 800—600 hPa. Compared to IASI, TES has a 
higher spectral resolution (0.06 vs 0.5 citT 1 ) and is thus more 
frequently sensitive to NH 3 in the boundary layer. Theoretical 
analysis indicates that the lower detection limit of the TES instru- 
ment for NH 3 is — 1 ppb (Shephard et al., 2011 ); it is likely 2—5 ppb 
for IASI (Clarisse et al., 2010). However, while TES has a smaller 
footprint (viewing area of approximately 5.3 km x 8.3 km), the 
observations are spaced far apart, which makes mapping difficult. 
On the other hand, IASI has much higher spatial coverage, 
providing global coverage twice daily with a 12 km x 12 km foot- 
print, which is important because NH 3 has high spatial variability. 
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Given the paucity of existing surface measurements of NH 3 , 
there is significant potential for enhancing our understanding of 
NH 3 using remote-sensing observations. It has been shown that TES 
is able to resolve spatial and temporal gradients in surface NH 3 
concentrations (Pinder et al„ 2011). Heald et al. (2012) used IASI 
measurements to study NH 3 concentrations over the U.S., finding 
underestimation of emissions in California and in the springtime in 
the Midwest. Zhu et al. (2013) have shown that TES observations of 
NH 3 can be used to constrain monthly average emissions. Model 
optimization for the U.S. improves the representation of measured 
NH 3 in April and October. Overall, remote-sensing constraints 
suggest that the NH 3 emission inventory is broadly underestimated 
in many parts of the U.S., particularly in the West. 

3.7. Particulate matter 

Aerosol data products from a subset of satellite sensors have 
been explored to estimate ground-level particle characteristics over 
land. MODIS and MISR measure particle abundance and composi- 
tion with nearly global coverage at moderate spatial resolutions, 
even over complex land surfaces. One of the most robust aerosol 
parameters that MISR and MODIS can retrieve is aerosol optical 
depth (AOD), also known as aerosol optical thickness (AOT). AOD is 
defined as the integral of aerosol extinction coefficients along the 
vertical atmospheric column from Earth’s surface to the top of the 
atmosphere. There are MODIS instruments aboard both the Terra 
(equator crossing time 10:30 am) and Aqua (equator crossing time 
1 :30 pm) satellites. A comparison of early- version Terra MODIS and 
Aqua MODIS AOD values indicated that differences vary by region 
without consistent patterns (Ichoku et al., 2005). 

The MODIS operational "Dense Dark Vegetation" (DDV) 
retrieval algorithm is designed to infer clear-sky aerosol properties 
over land surfaces that have low surface reflectance. This algo- 
rithm removes the surface impact by using the 2.1 pm band to 
estimate the surface reflectance, assuming a relationship of sur- 
face reflectance between this band and the blue (440 nm) and red 
(670 nm) bands as a function of both NDVI and scattering angle 
(Remer et al., 2005). AOD can then be derived by searching pre- 
computed lookup tables corresponding to different aerosol 
models. The MODIS operational algorithm has developed consid- 
erably over the last several years, with Collection 6 being recently 
described (Levy et al., 2013). 

The DDV algorithm is not performed over bright surfaces (i.e., 
when the 2.1 pm surface reflectance is above 0.15). The MODIS Deep 
Blue (DB) algorithm alleviates the bright-surface problem by 
employing radiance measurements from the blue channels to infer 
the properties of aerosols (Hsu et al., 2004). A maximum-likelihood 
method is used to match the appropriate values of AOD and SSA 
from a lookup table to the measured radiances. Studies over the 
Sahara Desert and nearby Arabian Peninsula, as well as the Gobi 
Desert, indicate that this algorithm is particularly sensitive to 
mineral dust over their bright source regions (Hsu et al., 2006; 
Fischer et al., 2009). Currently, MODIS DDV and DB AOD data 
(Collection 5.1 ) are reported as separate parameters in the Level 2 
data products (MOD04 and MYD04) at 10 km x 10 km nominal 
spatial resolution. 

MISR employs a unique multiangle design that allows it to 
observe the atmosphere through different effective path lengths, 
leading to a very different aerosol retrieval algorithm (Martonchik 
et al., 1998). To effectively separate surface-leaving light reflection 
from atmosphere-leaving light reflection over land, MISR uses the 
presence of spatial contrasts within a 17.6 km x 17.6 km retrieval 
region to derive an empirical orthogonal function (EOF) represen- 
tative of the region-averaged, surface-leaving light reflection. MISR 
defines a set of aerosol mixtures to represent aerosol types globally, 


which carry information about aerosol size distribution, shape, 
refractive index, and scale height. As a research instrument, MISR 
provides a rich aerosol dataset to study long-term spatial and 
temporal trends of particle mass, composition, and other infor- 
mation. MISR-retrieved AOD is reported in its Level 2 aerosol data 
product (MIL2ASAE) on a 17.6 km x 17.6 km grid. 

The stereoscopic nature of MISR data also enables it to retrieve 
the top heights of distinct plumes through a stereo-matching 
technique. Due to its narrow swath (~400 km), however, pre- 
liminary studies have shown that MISR only observes approxi- 
mately 10% of large wildfire smoke plumes (Kahn et al., 2008). 
Plume recognition is more likely to succeed close to the sources of 
fires, dust storms, volcanic eruptions, etc. The MISR plume-top 
height information is reported in the Level 2 cloud product 
(MIL2TCSP) on a 1.1 km x 1.1 km grid. 

MODIS and MISR AOD data provide the possibility for large- 
scale air pollution monitoring over land. Since the launch of 
MODIS and MISR, numerous studies have explored the relationship 
between satellite-retrieved aerosol parameters and ground-level 
air pollution measurements (e.g., Chu et al„ 2003; Liu et al„ 2005, 
2007; Wang and Christopher, 2003). Most studies have focused 
on the relationship between AOD and surface PM 2 . 5 , although a few 
have considered PMi 0 as well, since measurements of PM 2.5 outside 
North America and Western Europe are rare. The analytical 
methods have evolved from semi-quantitative descriptions of air 
pollution patterns, through correlation and simple linear regres- 
sion, to more sophisticated multivariate spatial and temporal 
models (van Donkelaar et al., 2010; Kloog et al., 2011; Lee et al., 
2011c; Liu et al., 2012a). 

Both PM and AOD can be observed using visible imagery from 
geostationary orbit, offering high spatial (30 min x 30 min) reso- 
lution. Examples of geostationary platforms with AOD products 
include the U.S. Geostationary Orbiting Earth Satellite (GOES) se- 
ries, the Japanese CMS and MTSAT, and the European SEVIRI. These 
AOD retrievals typically use only one or two visible channels 
(Prados et al., 2007). The geostationary AOD products offer higher 
observation frequency but have larger errors than MODIS or MISR 
(Paciorek et al., 2008). The Suomi-NPP VIIRS instrument was 
launched in October 2011 and has similar capabilities to MODIS but 
with higher resolution (6 km at nadir). The retrieval algorithm is 
developing rapidly and offers exciting possibilities to understand 
aerosol sources. 

3.8. Carbon dioxide 

Observing CO 2 concentrations from space is not difficult. In fact, 
many of the space-borne instruments have the capability to measure 
CO 2 columns, either by detecting the attenuation or emission of 
infrared radiation. The problem is to separate the C0 2 flux arising 
from natural sources from the flux due to human activity. Natural 
carbon fluxes are on the order of 300 GtC yr , while the anthro- 
pogenic contribution is only about 2% of this value, or ~6 GtC yr -1 . 
This presents a great challenge to the measurement instrumentation 
and the design of satellite platforms. Global measurements with an 
accuracy of ~ 1 ppm C0 2 are needed at high spatial and temporal 
resolution in order to be able to detect anthropogenic C0 2 in a useful 
way. Measurements of C0 2 concentrations using SCIAMACHY, AIRS, 
and TES have all been reported, as well as from the dedicated 
greenhouse-gas satellite GOSAT. It is generally held that satellite 
observations alone are at present unable to provide the required 
observational power for detecting point-source emissions, but that 
the coupling of satellite data with a ground-based sensor array and 
accurate meteorological data might be able to succeed (JASON, 2011 ). 

Specification of the necessary attributes of a space-borne global 
CO 2 observational platform that could detect anthropogenic 
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emissions began more than a decade ago. These were established 
with the design, development, and operation of the OCO and 
GOSAT platforms in mind. A number of inversion and precision 
analyses have been conducted to quantify how well CO 2 column 
measurements would be able to estimate surface sources and 
sinks of CO 2 (e.g., Chevallier et al., 2007, 2009; Miller et al., 2007; 
Baker et al., 2010; Pillai et al., 2010). Houweling et al. (2004) first 
confirmed that rather precise satellite measurements ( ~ 1— 2 ppm) 
would be needed to constrain sub-continental carbon fluxes. 
Lauvaux et al. (2009) focused on the effects of transport uncer- 
tainty on inversion modeling of sources and sinks. All these 
studies pointed out both the opportunities and the great chal- 
lenges ahead. 

Buchwitz et al. (2007) first reported direct observations of at- 
mospheric C0 2 concentrations between 2003 and 2005 from 
SCIAMACHY. Good agreement with CarbonTracker data was found, 
and the analysis was able to detect the C0 2 seasonal cycle and, for 
the first time, the year-to-year increase of CO 2 , which agreed with 
the CarbonTracker increase to within 1 ppm yr . Subsequent work 
compared a new SCIAMACHY column CO 2 retrieval with four 
TCCON sites and CarbonTracker data (Reuter et al., 2011). System- 
atic differences were in the range of 0.2— 0.8 ppm with a single- 
measurement precision of 2.5 ppm. The validation confirmed the 
potential of SCIAMACHY to provide information about column CO 2 
in regions with few in situ measurements. Observations of C0 2 from 
TES have been reported (Kulawik et al., 2010; Nassar et al., 2011). 
Kulawik et al. (2010) demonstrated good agreement with a variety 
of observational datasets, including CONTRAIL aircraft data, surface 
sites, AIRS data, and CarbonTracker. Nassar et al. (2011 ) inferred C0 2 
fluxes using inverse modeling of the satellite and surface mea- 
surements jointly with GEOS-Chem; the quantification of global 
fluxes from oceans and the terrestrial biosphere was in good 
agreement with other studies. 

The first results from GOSAT were reported in 2009 by the 
Japanese team (Yokota et al., 2009) using SWIR reflected light at 
1.6 pm. Pending further instrument calibration and validation, 
Yokota et al. (2009) were able to observe inter-hemisphere differ- 
ences in CO 2 concentrations, but with an underestimation of ab- 
solute concentration values. Subsequently, GOSAT retrievals of 
column-average CO 2 during its first year of operation (April 
2009— July 2010) were compared with TCCON and showed very 
close correlation (Butz et al., 2011). The seasonal CO 2 cycle and 
general source/sink patterns were evident. Further validation of 
GOSAT retrievals against TCCON data by Oshchepkov et al. (2012) 
and against TCCON and GEOS-Chem model calculations by Cogan 
et al. (2012) have also showed close correlation after removal of 
observations contaminated by cloud and aerosol. Global maps of 
CO 2 distributions from GOSAT have now been released as a Level 3 
product (Hammerling et al., 2012), generally in agreement with 
carbon modeling work, except for South America in July and 
August, which may point to inaccuracies in flux estimates due to 
the sparseness of in situ measurements there. 

A prime objective of space-based observations of C0 2 is to be 
able to detect emissions from individual point sources, which has 
yet to be fulfilled. Bovensmann et al. (2010) derived requirements 
for detecting the C0 2 emissions from a typical power plant. Mini- 
mum requirements of a spatial resolution of 2 km x 2 km, a wide 
swath width (~500 km), and a precision of 2 ppm could, they 
assert, be achieved with existing technology. The statistical un- 
certainty of the retrieval was found to be in the range of 12—36% of 
the emissions of a mid-sized power plant, with wind speed and 
aerosol concentration being the main interferences. The concept 
developed by Bovensmann et al. (2010), which they termed Car- 
bonSat, could detect and quantify C0 2 emissions from large point 
sources and would therefore have the potential to provide 


independent verification of emission reductions under a future 
global agreement. 

Kort et al. (2012) demonstrated the potential to quantify C0 2 
concentration enhancements due to human activity at the scale of 
megacities, using GOSAT observations over Los Angeles 
(3.2 ± 1.5 ppm) and Mumbai (2.4 ± 1.2 ppm). They estimate that 
changes as small as 0.7 ppm could be detected (in Los Angeles) from 
GOSAT observations at 95% confidence level. McKain et al. (2012) 
demonstrated the capability to determine changes in C0 2 emis- 
sions from an urban area (Salt Lake City) with a ground-based 
monitoring network and a CTM. Changes in emissions of ±15% 
could be detected at 95% confidence level. It was argued that in- 
tegrated C0 2 column measurements, of the kind that could be 
provided by space-borne instruments, would greatly enhance the 
detection of C0 2 emissions from source regions. Using SCIAMACHY 
data, Schneising et al. (2013) estimated an enhancement of 
1.1 ± 0.5 ppm for the U.S. east coast. 

With the limitations of current satellites to detect C0 2 emissions 
at the required level of detail and the failure of OCO in 2009, eyes 
turn to future satellite capability. The second version of the Orbiting 
Carbon Observatory (OCO-2) is presently under development with 
a scheduled launch date still to be determined, but possibly 2014/ 
2015; it would provide the capability needed to characterize 
sources and sinks of C0 2 at high resolution, when combined with 
ground-based measurements and related observations. A C0 2 
measurement platform in geosynchronous orbit would have many 
benefits, including the ability to revisit individual source locations 
on a time scale of hours rather than the 15—20 days of current low 
Earth orbit satellites (JASON, 2011). 

4. Methods for emissions estimation and examples of 
applications 

In this section, we first review analytical methods for constrain- 
ing emission estimates from remote sensing observations (4.1 ). This 
is followed by discussion of how more-detailed information about 
sources can be obtained by increasing the resolution of the satellite 
retrievals (4.2) and by the use of satellite data as a proxy for emis- 
sions trends (4.3). Then we discuss the potential to gain quantitative 
information about different kinds of sources from space-borne ob- 
servations by highlighting results from applications related to 
anthropogenic point sources (4.4), anthropogenic area sources (4.5), 
natural point sources (4.6), and natural area sources (4.7). We discuss 
synergy with field campaigns in a final sub-section (4.8). 

4.1. Inverse modeling techniques 

The most straightforward approach to developing “top-down” 
constraints on emissions directly attributes changes in observed 
concentrations to changes in emissions. This is commonly done for 
evaluating point sources (Sections 4.2 and 4.4) or for estimating 
long-term trends (Section 4.3). Such estimates can be based on 
observations alone. For assessment of area sources or point sources 
for which the signals are not distinct, because of transport, chem- 
istry, or the atmospheric lifetime of the species in question, more 
formal inverse methods are used. 

Inverse modeling seeks to formally develop constraints on 
emissions by combining three types of information: a set of ob- 
servations such as satellite measurements fy), an estimate of the 
distribution of sources of the species in question (x a , referred to as 
the a priori, prior, initial, or background), and a modeled simulation 
of the observed quantity using an estimate of the emissions (H(x)) — 
all with associated error estimates (see, e.g., Sandu and Chai, 2011). 
Many of the approaches to constraining emissions stem from Bayes 
theory, which describes how to determine the distribution of a set 
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of parameters (in this case, emissions) consistent with estimates of 
the remotely observed quantity that are closest to the distribution 
of these observations, while simultaneously accounting for prior 
information regarding reasonable values of these parameters. In 
general, solving for this probability distribution in full is computa- 
tionally prohibitive. In the following sections we review common 
methods for inverse modeling that approximate the Bayesian 
approach at varying levels of detail. The methods are distinguished 
by the ways in which the models are used, how the desired spatial 
and temporal scales of the analysis warrant various levels of 
complexity, and how errors in the inversion inputs are determined 
and propagated throughout the analysis. 

A common method for constraining emissions of short-lived 
species using remote sensing data is to average a set of satellite 
measurements over a particular location and derive top-down 
constraints on surface emissions using a mass balance approach. 
Model simulations based on an a priori emissions estimate are used 
to simulate atmospheric column or profile concentrations. The ratio 
of the simulated to the observed concentration, weighted according 
to estimates of the error in the observations and simulations, is 
then used to directly adjust the model emissions (e.g., Martin et al., 
2003). In this approach, the influences of atmospheric transport 
and chemistry on the relationship between emissions in one grid 
cell and concentrations in another grid cell are assumed to be 
negligible, presuming that the short chemical lifetime of the spe- 
cies translates to a smearing length scale less than the width of the 
model grid cell (see, e.g., Palmer et al., 2003). These assumptions 
are most suitable for short-lived species (e.g., NO*) and for rela- 
tively coarse (>1°) model simulations. However, even for short- 
lived species, the influence of emissions in a coarse air-quality 
model may extend beyond the column directly overhead (Turner 
et al., 2012). Kernels, or smoothing representations of the emis- 
sion field, have been used to account for the influence of adjacent 
cells (Toenges-Schuller et al., 2006; Boersma et al., 2008a). The 
procedure can be iterated to account for nonlinear relationships 
between emissions and column concentrations or profiles, or 
transport from neighboring grid cells (e.g., Zhao and Wang, 2009; 
Ghude et al., 2013). Differences in concentrations from multiple 
observations at different times can also be used to constrain the 
mass balance (Lin et al., 2010b). 

For longer-lived species or finer model resolutions, more explicit 
treatment of chemistry and transport becomes necessary. Formally, 
these approaches seek to minimize a cost function, J(x), which is a 
scalar measure of the misfit of the model ( H(x )) compared to ob- 
servations (3/), weighted by the observational error covariance 
matrix (R), plus the departure of the emissions from their a priori 
estimates (x b ), also weighted by the a priori error covariance 
matrix (B): 

JW = l/2(H(x)-y) T R- 1 (H(x)-y) + l/2(x-x b ) T B- 1 (x-x b ) 

(1) 

Approaches taken to finding the set of emissions that minimize 
this cost function depends upon several factors. The direct 
(or analytical) solution, which leads to the best linear unbiased 
estimate of emissions for linear systems, can be directly found as: 

x = x b + K(y - H(x)) (2) 

where K is the gain matrix, 

K = BH T (HBH T + R) 1 (3) 

and H is a Jacobian matrix that relates changes in emissions to 
changes in model estimates of observations determined from a 


(linearized) sensitivity analysis of the atmospheric chemistry and 
transport model H. This equation can be evaluated in a single step 
for all observations (direct or analytic inversion) or broken up into 
sequential steps (Kalman Filter) to update x, x b , and B, as a function 
of time. 

One of the main challenges with either approach is determining 
the Jacobian matrix. For a large number of emission values, or a 
large number of observations, this can simply become computa- 
tionally intractable. A common practice is to consider emissions 
aggregated at large spatial scales and then directly evaluate the 
Jacobian through successive forward model calculations (e.g., Heald 
et al., 2004; Stavrakou and Muller, 2006; Kopacz et al„ 2009). A 
drawback with this approach is the subjective lumping of emissions 
across broad regions, i.e., aggregation error. To define the Jacobian 
at finer scales, several approaches can be considered. More- 
efficient sensitivity techniques, such as direct decoupled or 
adjoint methods, may also be used to evaluate the Jacobian (e.g., 
Napelenok et al., 2008). Konovalov et al. (2006) define the rela- 
tionship between emissions and concentrations through Monte 
Carlo sampling of linear statistical relationships. Mijling and van 
der A (2012) introduced a Kalman-filter algorithm specifically 
designed for rapid inversion of short-lived atmospheric constitu- 
ents at a spatial scale of ~ 25 km x 25 km, by approximating the 
Jacobian matrix through 2D trajectory analysis. 

An alternative approach is the ensemble Kalman filter method 
(Evensen, 1994). This does not require specification of Jacobian 
matrices; instead, the influence of emissions and their un- 
certainties on the model estimates is assessed through analysis of 
the spread of ensembles (typically 20—40) of forward model eval- 
uations. For example, Miyazaki et al. (2012a) use the ensemble 
Kalman filter approach to constrain NO* emissions with OM1 ob- 
servations. Such approaches, becoming increasingly common in 
weather forecasting systems, have become attractive in recent 
years given their scalability on parallel computing systems. 

Lastly, the minimum of the cost function may be sought itera- 
tively using gradient-based optimization techniques. This method, 
referred to in the context of air quality forecasting as 4D-Var, is 
based on the calculus of variation, wherein a perturbation in a 
scalar, receptor-based metric is propagated backward in time 
through a set of equations that are auxiliary (adjoint) to the locally 
linearized system. This affords calculation of the gradient with 
respect to numerous parameters simultaneously. These gradients 
are then used in conjunction with an optimization scheme to 
minimize the cost equation. The 4D-Var method was first used to 
constrain emissions in an Eulerian air-quality model of chemically 
active species in the troposphere by Elbern et al. (2000). Effectively, 
this method affords constraints on model emissions at the native 
model resolution (e.g., Kopacz et al., 2009; Jiang et al., 2011; 
Hooghiemstra et al., 2012), which minimizes the aggregation er- 
ror incurred when solving Equation (2) directly for a vector or 
emissions defined across broad areas. The approach does, however, 
require the existence of an adjoint model, and estimating the un- 
certainty in the resulting emissions requires additional effort. 

Within the frameworks for developing top-down emissions 
constraints described above, there are several challenges related to 
the use of new remote sensing constraints, applications to new 
inverse modeling objectives, and ascertaining the uncertainties in 
the resulting inventories. One issue at the forefront of inverse 
modeling is to use measurements from multiple species, or from 
multiple platforms. Data assimilation techniques have typically 
focused on the recovery of emissions of a single species. For 
example, Dubovik et al. (2008) used MODIS AOD to estimate fine- 
and coarse-mode aerosol emissions. More recently, Huneeus et al. 
(2012) demonstrated the simultaneous estimate of global emis- 
sions of multiple gaseous and aerosol species (dust, sea salt, BC, OC, 
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and S0 2 ) using a 4D-Var approach assimilating MODIS total and 
fine-mode AOD. Miyazaki et al. (2012b) have used the same data 
assimilation system to combine observations of four interrelated 
species — N0 2 , O 3 , CO, and HNO 3 — from four satellites— OM1, TES, 
MOP1TT, and MLS — with additional aircraft and ozonesonde data. 
More work is needed on top down-estimates for multiple species 
(such as N0 2 and VOC or N0 2 and O 3 ). Early attempts (Zhang et al., 
2008; Carmichael et al., 2008) produced mixed results, where in 
some regions the assimilation was pulled in opposite directions by 
the different data sets. 

There are also issues with trying to assimilate a single parameter 
from multiple retrievals. For example, Jones et al. (2009) compared 
an analytic inversion, Equation (2), of MOPITT and TES in the tropics 
in November 2004. While, in general, the a posteriori emissions 
showed some consistency, over North America — where in boreal 
fall emissions make only a small contribution to total CO abun- 
dance — MOPITT and TES adjusted initial estimates in opposite di- 
rections. Fortems-Cheiney et al. (2009) presented a first 
comparison of inversions using IASI and MOPITT data and found 
that the regional totals were broadly consistent. George et al. 
(2009) compared IASI CO retrievals with MOPITT, AIRS, and TES 
(see Fig. 3) and found, on average, total column discrepancies 
of ~ 7% in the Northern Hemisphere and in the equatorial region. 
All the studies concluded that there was a need for better charac- 
terization of biases among data sets. 

Model errors also pose a significant limitation on the accuracy of 
top-down emission estimates. Arellano and Hess (2006), in an 


IASI 



inter-model comparison, showed that the results of inversions are 
particularly sensitive to the model treatment of convection and to 
transport out of the boundary layer. A linearization error is intro- 
duced when the chemical loss of CO is accounted for by specified 
OH abundances (Muller and Stavrakou, 2005), an approach 
commonly used to minimize the computation cost of the inversion. 
Differences between inversions based on remote sensing observa- 
tions of column concentrations compared to surface observations 
can also be used to identify components of model transport error 
that impact the emissions constraints (Jiang et al., 2013). 

The selection of source types to be constrained also impacts the 
inversion solution. For example, since CO observations mainly 
constrain total emissions, an increase in one source category may 
be compensated for by a decrease in another. A careful character- 
ization of the state vector using the spatiotemporal resolution and 
precision of observations rather than geographical boundaries is 
required. Accounting for the CO contribution that is not directly 
emitted adds another level of uncertainty. In most studies (e.g., 
Palmer et al., 2003; Heald et al., 2004; Arellano and Hess, 2006; 
Jones et al., 2009; Kopacz et al., 2009), the biogenic NMVOC 
source, representing a main source for CO precursors, is aggregated 
into a background global chemical source of CO. Stavrakou and 
Muller (2006) and Jiang et al. (2011) explicitly solved for the 
biogenic NMVOC source at the same spatial resolution at which 
the combustion CO emissions were constrained and showed 
that aggregating the biogenic NMVOC source into the global 
background could result in significant over-adjustment of the a 
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Fig. 3. Averaged IASI, MOPITT, AIRS, and TES CO total column distributions, binned on a 1° x 1° grid for August 2008. All observations are for daytime and were cloud-filtered 
following the recommendations provided by each retrieval team (adapted from George et al., 2009). 
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posteriori CO emissions. The range of solutions from different in- 
versions was shown to be much larger than the a posteriori errors, 
indicating that the solution is more sensitive to modifications to the 
inverse system than a posteriori errors suggest. Without a reason- 
able model error estimate, the results of inverse analyses will be 
sensitive to the choice of the observation domain and the proximity 
of observations to the source location of interest. 

4.2. Oversampling/spatial smoothing 

Some species have well-defined plumes that tend to be 
smoothed out during spatial averaging. It is possible to improve the 
definition of these plumes by oversampling the original swath data 
to a higher resolution. Oversampling is a powerful tool, because 
higher resolution data are produced from coarse observations and 
noise is reduced, thereby enhancing weak signals. This technique 
has been widely used for NO 2 (e.g., Martin et al., 2003; Beirle et al., 
2004b; Lamsal et al., 2008) and was recently applied to high-quality 
S0 2 pixels from the OM1 sensor to a resolution of 3 km x 3 km 
for retrievals above Mexico City (de Foy et al., 2009). Using this 
technique, the separate plumes and dispersion patterns of an in- 
dustrial complex and a passively degassing volcano became clearly 
visible. By averaging OMI N0 2 columns to a grid resolution of 
0.025° x 0.025° over California, Russell et al. (2010) were able to 
identify spatial patterns in the weekday-weekend concentrations 
in the Los Angeles basin. 

Fioletov et al. (2011) used a similar technique to evaluate the 
strengths of large point sources of SO 2 in the U.S. on a 2 km x 2 km 
grid. Because the elevated S0 2 values associated with a point source 
generally become insignificant beyond ~ 50 km, they averaged all 
OMI pixels within a 12 km radius of the source. At the highest 
spatial resolution currently available from OMI (>300 km 2 ) the SO 2 
detection limit is 0.1— 0.3 DU, which corresponds to annual emis- 
sions larger than ~70 Gg SO 2 yr . Currently, there are about 40 
point sources in North America (smelters and coal-fired power 
plants in Canada, the U.S., and Mexico) that exceed this source 
strength. Fioletov et al. (2011) also observed reductions of column 
SO 2 of ~40% that were consistent with reported decreases in 
emissions between 2005—2007 and 2008—2010. Fig. 4 shows the 


high-resolution S0 2 retrievals reported in Fioletov et al. (2011). 
Applying the same method to Alberta, Canada, McLinden et al. 
(2012) showed annual increases in N0 2 columns between 2005 
and 2010 due to the development of oil sands (see Section 4.5). 
These techniques have generated higher-resolution emissions dis- 
tributions that significantly improve our ability to detect and 
quantify emission sources. 

4.3. Timely updates to emission trends and other temporal 
variations 

The bottom-up approach of compiling emission inventories by 
aggregating activity data and emission factors provides the most 
comprehensive treatment of emissions, but this approach is a major 
undertaking that can require time to complete. For example, at the 
time of writing (2013) the current NEI from the U.S. EPA is for the 
year 2008, though estimates for more recent years are available in 
other products (see below). More timely information could 
improve the accuracy of model simulations for time periods after 
the release of the inventory. Satellite observations offer near-real- 
time information on a range of species, as discussed above, and 
the trends in the satellite measurements can be used to update 
bottom-up emission trends, pending the gathering of source in- 
formation to construct a new bottom-up inventory. 

Lamsal et al. (2011) demonstrated that satellite observations of 
NO 2 provide timely information about changes in anthropogenic 
NO* emissions, when combined with model information on the 
relationship between changes in tropospheric NO 2 columns and 
changes in NO* emissions. Fig. 5 shows the spatial variation of 
anthropogenic NO* emissions for 2005 from a bottom-up inventory 
and a predicted NO* inventory for 2010 from OMI measurements. 
The predicted NO* inventory (24.8 Tg N yr -1 ) is 2.5% higher than the 
bottom-up inventory (24.2 Tg N yr -1 ) for regions dominated by 
anthropogenic NO* emissions. The difference between the pre- 
dicted inventory for 2010 and the bottom-up inventory for 2005 
presented in Fig. 5 (bottom frame) shows significant reductions 
over the eastern U.S. .Japan, and parts of Europe, and increases over 
eastern China. Changes in anthropogenic NO* emissions during 
2005—2010 indicate a decrease of 22.7% and 13.7% over the U.S. and 
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Fig. 4. Mean S0 2 burdens over the Ohio River Basin for 2005-2007 (left) and 2008-2010 (right) measured by OMI, confirming a substantial reduction in S0 2 pollution around the 
largest coal-fired power plants, as a result of the implementation of S0 2 emission control measures (adapted from NASA Earth Observatory, as reported in Fioletov et al., 2011). 
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Fig. 5. Spatial distributions of anthropogenic NO, emissions at 1° x 1.25°: bottom-up inventory for the year 2005 (top); inventoiy predicted from OM1 N0 2 observations for the year 
2010 (middle). Inventory totals refer to colored regions only. The difference between the OMI-derived emissions for 2010 minus the bottom-up emissions for 2005 is shown at the 
bottom (adapted from Lamsal et al., 2011 ). 


OECD Europe, respectively, and an increase of 24.9% in East Asia, 
with a 26.8% growth in Chinese NO* emissions. The predicted 
decrease of 22.7% for the U.S. is very similar to the decrease 
calculated from U.S. EPA trends data for 2005 and the unreleased 
2010 NEI, providing confirmation of the updating approach. 

A similar approach may be applicable to other species such as 
CO and NH3 to monitor emission changes over time (e.g., Worden 


et al., 2013). Fortems-Cheiney et al. (2011) conducted an inversion 
of the 10-year time series of MOP1TT CO and interpreted interan- 
nual changes in terms of variability due to biomass burning, 
climate, and socioeconomic factors. However, this study used 
MOPITT v4 data, which did not account for instrument drift that has 
resulted in a 0.8 ± 0.2% yr _1 increase in the retrieved column 
(Deeter et al., 2010). Recent MOPITT v5 data (Deeter et al., 2012) 
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account for this drift and compare well with temporal trends from 
other satellite data, as shown in Fig. 6, which is from a study by 
Worden et al. (2013). This study did not analyze changes in emis- 
sions directly, but compared trends in retrieved CO columns from 
MOPITT, AIRS, TES and IASI for 2000-2011. They found a modest 
decreasing trend ( — 1% yr _1 ) in total column CO over the North- 
ern Hemisphere and a less significant, but still decreasing trend in 
the Southern Hemisphere. The decreasing trend over Europe and 
the U.S. has been reported in other studies and is related to re- 
ductions in anthropogenic emissions (Zellweger et al., 2009; 
Angelbratt et al., 2011; EPA, 2013b). 

Not only can time series of satellite retrievals help to extend and 
interpret long-term trends in emissions in various parts of the 
world (e.g., van der A et al., 2008; Lin et al., 2010a; Schneider and 
van der A, 2012), but the effects of short-term changes in emis- 
sions (days to weeks to a year) — due to deliberate human inter- 
vention or to autonomous forces — can be observed. This suggests 
significant capability to improve the sensitivity of emission in- 
ventories by convolving satellite-derived temporal variability. 

Anomalies in annual trends have been reported in a number of 
studies, both in the U.S. and other parts of the world. Russell et al. 
(2012) noted from examination of NO 2 retrievals in California that 
there had been a decrease in diesel truck activity attributed to the 
economic recession (Russell et al„ 2012). Castellanos and Boersma 
(2012) observed similar N0 2 reductions over Europe during 


2008—2009, though in both California and western Europe the 
continued implementation of NO* emission controls makes it 
difficult to assign causation unequivocally. Satellite-detected effects 
of the economic recession on NO* emissions have also been re- 
ported for China (Lin and McElroy, 2011), Greece (Vrekoussis et al., 
2013), and international shipping (de Ruyter de Wildt et al., 2012). 
And in a related shorter-term study of economic impact, Yoshida 
et al. (2010) showed that the effects of Hurricanes Katrina and 
Rita in 2005 could be detected in the form of reduced OMI N0 2 
signals around the Gulf of Mexico, associated with disruptions to oil 
and gas production, oil refining, and power generation. 

Changes in emission levels on the scale of days to months due to 
particular policy interventions or other driving forces have been 
observed. The first reported example of this was the change in NO* 
emissions in Beijing at the time of the Sino-African Summit in 
November 2006 (Wang et al„ 2007). Beijing municipal authorities 
instituted policy measures to limit traffic in Beijing during a six-day 
period, with an estimated 30% reduction in the numbers of vehicles 
on the road. OMI observed a reduction in column densities of about 
a factor of two during the period of the Summit. In some ways, this 
event was a dress rehearsal for the 2008 Beijing Olympic Games, in 
anticipation of which even more drastic emission reductions were 
implemented, covering both mobile and stationary sources, and 
extending from Beijing itself to several surrounding provinces. This 
unprecedented intervention to reduce emissions over roughly a 
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Fig. 6. Comparison of MOPITT, AIRS, TES, and IASI retrieved annual average CO column amounts for the period 2002-2012 over four different world regions: E. China, E. USA, 
Europe, and India (adapted from Worden et al., 2013). 
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two-month period stimulated many research studies utilizing 
ground-based measurements, satellite observations, and emission 
inventories: N0 2 from GOME-2 and OMI (Mijling et al., 2009), N0 2 
and S0 2 from OMI and CO from MOPITT (Witte et al., 2009), aero- 
sols from MODIS (Cermak and Knutti, 2009; Liu et al., 2012a), 
aerosols from MODIS and N0 2 from OMI (Zhang et al., 2009b), and 
CO from MOPITT (Worden et al., 2012a). In all cases, the satellite 
measurements confirmed the emission reductions that were re- 
ported to have been implemented. A return to pre-existing emis- 
sion levels after the conclusion of the Olympic Games was in some 
cases masked by the onset of an economic recession in the fall of 
2008. The experience of the Beijing Olympic Games was revisited 
during the 2010 World Expo in Shanghai, when Hao et al. (2011) 
reported significant (8—14%) emission reductions of N0 2 , CO, and 
aerosols during the six-month period of the Expo. 

Because of the daily coverage provided by some of the satellites, 
it is possible to investigate even finer temporal cycles of emissions. 
In the first study of its kind, Beirle et al. (2003) were able to discern 
different weekly cycles of NO* emissions in different parts of the 
world from GOME N0 2 columns. In countries with a largely 
Christian tradition, there was a distinct reduction in emissions on 
Sundays and a slight reduction on Saturdays. In Israel the minimum 
occurred on Saturdays, and in Islamic cities of the Middle East there 
was a slight reduction on Fridays. In China, no weekly effect was 
observed, possibly due to the fact that emissions are dominated by 
power plants and industrial plants that operate seven days a week. 
Kaynak et al. (2009) extended the work of Beirle et al. (2003) by 
using SCIAMACEIY retrievals during the period 2003—2005 over 
rural, urban, and power-plant areas of the continental U.S. A 
distinct weekly pattern was observed in cities, with a minimum on 
Sundays and a smaller reduction on Saturdays, as found by Beirle 
et al. (2003). Rural areas showed no clear weekly pattern. Areas 
dominated by power plants showed only a small Sunday reduction. 
It was possible to infer the relative contributions of mobile sources 
to total emissions from these weekly patterns, with results differing 
somewhat from what would be obtained from bottom-up in- 
ventories, suggesting a possible application of the satellite data to 
refine day-to-day variability of emissions over urban areas. Russell 
et al. (2010) examined the situation over California from OMI data 
and confirmed that weekday emissions are higher than weekend 
emissions. They observed that transportation studies have ascribed 
much of the difference to the absence of heavy-duty diesel trucks 
on weekends. Similar weekly cycles were reported over Europe by 
Zhou et al. (2012). Even observations of AOD have revealed weekly 
cycles over the U.S., central Europe, India, and the Middle East (Xia 
et al., 2008). 

Since each NO* source has specific characteristics that deter- 
mine the seasonal variation of the tropospheric N0 2 columns, the 
seasonality of N0 2 satellite observations can be used to identify the 
dominant sources of emissions (van der A et al., 2006, 2008; Ghude 
et al., 2010; Zhang et al„ 2012). The idea builds on seasonal char- 
acteristics that have been exploited in earlier studies to understand 
sources. For example, regions that are dominated by anthropogenic 
emissions have a seasonal N0 2 maximum in the winter (Martin 
et al., 2003); biomass burning of forests and savanna usually 
takes place during the dry season such as early spring; and soil 
emissions can result in enhanced N0 2 signals in summertime 
(Jaegle et al., 2004). Using these characteristics, van der A et al. 
(2008) first identified the dominant sources of NO* emissions on 
a global scale at a resolution of 1° x 1° on the basis of GOME and 
SCIAMACHY measurements. On a regional scale, they found that 
eastern and western China are dominated by anthropogenic and 
natural emissions, respectively. Recently, Zhang et al. (2012) 
applied a similar methodology to a multiyear analysis from 1996 
to 2010. They observed that the highly polluted regions (areas 


dominated by anthropogenic emissions) in China have expanded 
from the east to the central and the west, and new highly polluted 
regions have formed throughout the nation in the past 15 years. 
This technique shows the potential for determining changes over 
time in the spatial extent of different area sources having different 
emitting patterns that can then be fed back into bottom-up 
inventories. 

4.4. Anthropogenic point sources 

Emissions from large, man-made point sources are clearly one of 
the major targets for measurement from space, both because they 
represent a strong signal that should be within the detection limits 
of available instruments and also because they are potentially large 
polluters that command the attention of environmental regulatory 
bodies. The ability to detect changing emissions from these types of 
sources has the potential to aid in the verification of region-wide 
pollution control policies and in determining the compliance of 
individual point sources with emission control requirements. 
Consequently, anthropogenic point sources, particularly of NO* and 
S0 2 , have been the subject of many research studies. For application 
to U.S. air quality management, the interplay between observed 
emissions and regulatory policy is a key component of such work. 

During the era of measuring air quality from space ( ~ mid- 
1990s to present), there have been several important environ- 
mental efforts in the U.S. directed at reducing the emissions of NO* 
and S0 2 from coal-burning power plants. For example, EPA issued 
the 1998 NO* State Implementation Plan (SIP) Call with the intent 
to reduce NO* emissions during the summer season and thereby to 
reduce ground-level O 3 concentrations. And in 2005 EPA issued the 
Clean Air Interstate Rule (CAIR) with the goal of decreasing NO* and 
S0 2 emissions even further. Both the SIP Call and CAIR apply to the 
eastern U.S., where emissions in an upwind state can significantly 
impact air quality in another, downwind state. 

There are a number of methods for significantly reducing 
emissions of NO* and S0 2 from the effluent of stationary sources, 
most importantly Selective Catalytic Reduction (SCR) for NO* and 
flue-gas desulfurization (FGD) for S0 2 . SCR and FGD units can 
effectively remove up to 90% of NO* and S0 2 , respectively, from the 
effluent stream. Application of these systems on a number of power 
plants, coupled with other strategies to limit emissions from 
polluting plants, has caused NO* and S0 2 emissions from power 
plants to decrease substantially from 1999, when about one- 
quarter of total NO* emissions and almost all S0 2 emissions in 
the U.S. were associated with power generation. Given the short 
lifetimes of NO* and S0 2 , data collected from space can be used to 
confirm reductions reported by individual power plants, particu- 
larly those not impacted by large area sources such as cities. 
Inferring emissions or changes in emissions from point sources 
using satellite data requires careful attention, however, as the es- 
timate depends on the assumed lifetimes for S0 2 and N0 2 , which 
can be highly variable due to weather and other factors (see, e.g., 
Beirle et al., 2011; Walter et al., 2012; Zhou et al., 2012). 

The first work on U.S. power plants was reported by Kim et al. 
(2006), who used GOME and SCIAMACHY data to examine the 
long-term evolution of N0 2 columns in the eastern U.S. from 1997 
to 2005, with the intention of assessing regional changes in NO* 
emissions associated with the 1998 NO* SIP Call. They inferred that 
NO* emissions from power plants in the Ohio River Valley had 
decreased by about 35% over this period, which was in line with 
reported measured emission changes from Continuous Emissions 
Monitoring Systems (CEMS). This result was in contrast to the 
Northeast urban corridor, where there was no clear trend in 
emissions. In subsequent work Kim et al. (2009) reported an 
analysis of 2005 N0 2 satellite data for 13 individual power plants in 
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the western U.S. Agreement with WRF-Chem modeled columns 
was found to be excellent for three retrievals: SCIAMACHY 
(r = 0.93); the University of Bremen OMI (r = 0.97), and NASA OMI 
(r = 0.95). These plants are considerably more remote from other 
sources than plants in the eastern U.S., making the signals from 
their plumes more distinctive. Russell et al. (2012) used NO 2 col- 
umn data from the OMI BEHR retrieval algorithm to infer that NO* 
emission changes from 23 large power plants throughout the U.S. 
were variable because of region-specific regulations, decreasing on 
average by 26 ± 12% from 2005 to 2011. 

China and India are excellent test beds for the study of point 
sources. Both countries contain large coal-fired power plants that 
are relatively poorly controlled and therefore have high emission 
rates. Not only that, the boom period of the mid-2000s saw the 
construction of many new plants of this kind, meaning that we 
have the luxury of being able to observe the satellite responses 
before and after the sources of emissions came into existence. This 
is rarely possible in the developed world these days. The close 
relationship between N0 2 columns and the locations of large coal- 
fired power plants in China and India is shown in Fig. 7. 

The potential for investigating emissions from power plants in 
China was triggered by the groundbreaking paper of Richter et al. 
(2005), which reported the large growth (~50%) in NO* emis- 
sions in China inferred from GOME and SCIAMACF1Y measurements 
during the period 1996—2004. This was confirmed by a combined 


bottom-up and top-down examination of China’s NO* emission 
trends over the same period (Zhang et al., 2007). Subsequently, 
Zhang et al. (2009a) studied the NO x emissions from new power 
plants built in the period 2005—2007 in Inner Mongolia, China, 
where the power plants are relatively far removed from large urban 
areas. OMI summertime N0 2 columns were found to be closely 
correlated with emissions from the new large sources of pollution. 
Quantitative agreement (rates of increase) were very good in cases 
where new sources were added to already-developed regions 
(R 2 = 0.94), but not so good in remote regions, due to uncertainty in 
the background NO* ( R 2 = 0.73 for all plant locations). This work 
was subsequently expanded to a wider set of Chinese power plants 
with similar results (Wang et al., 2010, 2012). 

Li et al. (2010b) extended the study of Chinese power plants to 
S0 2 , with the aim of demonstrating that OMI could detect the 
presence or absence of S0 2 emission control equipment and 
possibly even quantify its efficiency. Following Zhang et al. (2009a), 
the S0 2 /N0 2 column ratios were determined for a number of power 
plants in Inner Mongolia in 2005—2008, during which period the 
plants were scheduled to install FGD equipment. Dramatic declines 
in S0 2 were observed at some plants — implying the implementa- 
tion of pollution controls — while the N0 2 signal remained rela- 
tively constant — implying steady electricity generation over the 
time period. Li et al. (2010b) clearly demonstrated that there is 
great potential to verify the effectiveness of pollution measures 
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Fig. 7. Locations of thermal power plants in Asia, particularly in China and India, in relation to 2005-2010 average OMI N0 2 column density measurements. Note that the OMI scales 
are different in the two parts of this image. 
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from space, which was bolstered by the high level of interest in this 
work expressed by China’s Ministry of Environmental Protection. 

There have been fewer papers devoted to power plants in India, 
though the situation is similar to China in that many new coal- 
fired plants with relatively poor controls were built after the 
availability of high-performance satellite measurements. In many 
ways India is a more desirable target than China because of shorter 
NO 2 lifetimes in hotter conditions, smaller solar zenith angle, 
weaker seasonality, and lower area emissions around power 
plants (Lu and Streets, 2012). GOME and SCIAMACHY trends over 
India for 1996—2006 have been used to identify major hotspots of 
NO* emissions, some of which were power plants (Ghude et al., 
2008). Prasad et al. (2012) used OMI data to study power plants 
in India over the period 2005—2007. Lu and Streets (2012) com- 
bined results from four satellite instruments, GOME, SCIAMACHY, 
OMI, and GOME-2 to quantify the growth in NO„ emissions from 
Indian power plants. 

Overall, it can be concluded that emissions of NO* and SO 2 from 
coal-fired power plants represent promising targets for application 
of satellite observations. The China and India experience needs to 
be transferred to the North American context with detailed, unit- 
specific comparison of space-based emission estimates with 
CEMS measurements under a variety of location, seasonal, and 
meteorological conditions. An inter-comparison of derived re- 
lationships between, for example, power-plant NO* emission rates 
and NO 2 column concentrations from different studies in different 
parts of the world would also be a step forward. 

Metal smelters are the largest anthropogenic point sources of 
S0 2 and have therefore been targets of study to test the capability of 
space-borne instruments. Carn et al. (2007) first reported the 
detection of S0 2 emissions from two copper smelters in Peru using 
OMI. Annual SO 2 emissions from the larger of the two smelters, Ilo, 
were estimated at 0.3 (0.2— 0.5) Tg. Interestingly, the emission rate 
appeared to decrease by about 40% between late 2004 and early 
2005, which the authors of the paper speculate may have been 
caused by modernization of the plant and an increase in the SO 2 
capture rate. This mirrors the work on S0 2 capture by power plants 
in China (Li et al., 2010b). Khokhar et al. (2008) analyzed SO 2 
emissions from the same two smelters in Peru and the large smelter 
complex at Norilsk in Russia, using a seven-year time series of 
GOME data. They estimated annual S0 2 emissions from the Ilo 
smelter at 1.095 ± 0.2 Tg, larger than the value obtained by Carn 
et al. (2007). Uncertainties related to the assumed S0 2 lifetime 
and AMF are cited as possible reasons for discrepancies. 

The copper and nickel smelters at Norilsk in Russia are 
acknowledged to be the largest anthropogenic point sources of SO 2 
in the world. Khokhar et al. (2008) estimated their annual emis- 
sions to be 1.685 ± 0.3 Tg from GOME data. Walter et al. (2012) were 
able to reassess Norilsk emissions using the DOAS system on board 
the CARIB1C aircraft, yielding an estimate of annual SO 2 emissions 
of 0.92 ± 0.5 Tg. The estimate derived from closely synchronized 
OMI measurements was 0.7 (0.6— 0.9) Tg. Walter et al. (2012) report 
on several bottom-up estimates from such sources as the EDGAR 
inventory and the Norilsk Nickel Company itself, which suggest 
that emissions “ought” to be higher than the satellite observed 
amounts — in the region of 2 Tg yr '. 

In a similar study Carn et al. (2004) observed the S0 2 that was 
released from a month-long fire at the Al-Mishraq State Sulfur 
plant in Iraq in 2003. Using TOMS data they estimated a release 
of ~0.6 Tg of S0 2 in total — roughly consistent with the reported 
inventory loss of sulfur— which made it the largest non-volcanic 
SO 2 emission event observed to date. Kearney et al. (2009) subse- 
quently resampled the MODIS data to the spatial resolution of 
TOMS and reanalyzed the plume S0 2 for a single day, 29 June 2003, 
finding good agreement between the two data sets. Large 


anthropogenic S0 2 point sources like these are clearly visible in the 
satellite retrievals, but, fortunately, they are few and far between. 

4.5. Anthropogenic area sources 

Anthropogenic area sources are a natural target of opportunity 
for space-based emissions estimation, as they are individually small 
in size, spatially dispersed, and difficult to accurately represent in 
emission inventories. Unconventional energy extraction processes 
are such an area source. The important work of McLinden et al. 
(2012) examined emissions of NO* and S0 2 from oil/tar sands op- 
erations in the Canadian province of Alberta, using the high- 
resolution (~30 km x 50 km) retrievals previously discussed in 
the context of Fioletov et al. (2011) and Boersma et al. (2011). In- 
creases in N0 2 columns between 2005 and 2010 were estimated to 
be a remarkable 10.4 ± 3.5% yr , reflecting the dramatic growth in 
exploitation of this resource both in intensity and spatial extent. 
Similar applications to such activities as oil and gas extraction in the 
Colorado Northern Front Range and other locations are likely to 
follow. 

Shipping emissions represent an area source category that has 
posed problems for bottom-up inventories due to the wide variety 
of source types and poor statistics on numbers of vessels and 
emission factors. The first remote detection of ship tracks was re- 
ported by Beirle et al. (2004a) along the international shipping lane 
from Sri Lanka to the Strait of Malacca across the Indian Ocean. They 
used a 6-year (1996— 2001) composite of GOME N0 2 data. Seasonal 
differences in wind direction were used to derive a mean lifetime of 
NO* in ship exhaust of 3.7 (1.9— 6.0) hours. Using this result, Beirle 
et al. (2004a) derived an estimate for annual shipping emissions in 
the Indian Ocean domain of 23 (10—73) Gg N, in good agreement 
with inventory estimates. In that same year, Richter et al. (2004) 
reported SCIAMACHY measurements of N0 2 associated with ship- 
ping not only in the Indian Ocean, but also in the Red Sea and the 
South China Sea. Reasonable agreement with shipping emission 
inventories was obtained, though in some locations uncertainty 
was high. Better knowledge of the N0 2 lifetime was needed — they 
used a value of 5.6 h, longer than Beirle et al. (2004a) — in order to 
more accurately characterize the horizontal and vertical dispersion 
of the plumes and NO* removal. 

This work was later extended by Franke et al. (2009) using both 
SCIAMACHY and GOME-2 to obtain monthly information about 
shipping emissions, though springtime biomass burning over India 
and autumn burning over Indonesia interfered with the analysis. 
They concluded that an annual emission estimate of 90 Gg N was 
reasonable, increasing at a rapid rate that could account for dif- 
ferences between GOME and SCIAMACHY measurements over the 
period 1996—2007. Marbach et al. (2009) detected HCHO from 
shipping in the Indian Ocean using GOME data. The HCHO is not 
emitted directly from the ships but is thought to be formed by some 
unknown photochemical mechanism in or close to the ship plumes. 
Most recently, de Ruyter de Wildt et al. (2012) used data from four 
space-borne instruments to observe N0 2 columns from shipping 
activity in a region stretching from the Mediterranean Sea to the 
South China Sea. SCIAMACHY detected four major shipping lanes: 
the Mediterranean Sea between Italy and Tunisia, the Red Sea, the 
Indian Ocean (as reported by Beirle et al., 2004a), and the South 
China Sea leading northeast from Singapore to Chinese ports. They 
developed time trends in N0 2 column density from 1996 to 2010 
and compared them with trends in shipping cargo volume and 
international trade volume. 

One major potential application for satellite observations is to 
estimate the emissions in a large metropolitan area, where there 
may be very many small sources of different types that in aggregate 
represent a large emissions contribution spread over a wide 
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domain. In the U.S., N0 2 urban signatures and their weekly cycles 
were first reported for five cities (New York, Los Angeles, Chicago, 
Dallas, and Detroit) using GOME measurements (Beirle et al., 2003). 
Later work addressed 14 western cities (Kim et al., 2009), including 
annual trends and weekly cycles. Similar work was then performed 
for several cities located across the U.S. (Kaynak et al., 2009) and 
cities in California (Russell et al., 2010). Satellite observations show 
that N0 2 columns in U.S. cities have generally been declining during 
the observational record, which is consistent with the imple- 
mentation of the requirements of the Tier 1 and 2 vehicle standards 
and the gradual turnover of the fleet of light-duty vehicles. Russell 
et al. (2012) used the BEHR OMI N0 2 retrieval data to estimate an 
average total reduction of 32 ± 7% for U.S. cities from 2005 to 2011 
with a 34% decrease in N0 2 from mobile sources. U.S. EPA trends 
data show a reduction of 31% between 2005 and 2011. Russell et al. 
(2012) attributed part of the observed decline in OMI N0 2 columns 
to the turnover in the mobile source fleet and part to the economic 
recession that began in 2008. Specifically, they reported that N0 2 
columns decreased by 6 ± 5% yr _1 before the recession and 8 ± 5% 
yr _1 during the recession. 

Developing a bottom-up emission inventory for a megacity in 
the developing world is an even more daunting task. Conventional 
techniques for doing so are confounded by influxes of migrant 
workers and unconventional vehicles from the countryside, never- 
ending construction activities with their ubiquitous diesel en- 
gines, unregistered vehicles, small industrial operations, household 
burning of solid fuels, etc. These difficulties can, to a certain extent, 
be by-passed by a satellite observation, which sees only the inte- 
grated product of multifarious small sources. Indeed, the potential 
of satellite N0 2 observations to yield information about the size, 
shape, land use, and dynamics of an urban area has prompted a new 
breed of demographic/socioeconomic studies, exemplified by 
Bechle et al. (2011) and Novotny et al. (2011). For urban studies 
outside North America, the reader can consult satellite analyses of 
cities worldwide (Clerbaux et al., 2008; Kar et al., 2010), in Europe 
and the Middle East (Konovalov et al., 2010), in the eastern Medi- 
terranean (Kanakidou et al., 2011), and in southeastern Europe 
(Zyrichidou et al., 2009). Beirle et al. (2011 ) used OMI measurements 
over eight large urban agglomerations to probe the lifetime of NO„, 
with an eye to elucidating the formation of urban O3. A good rela- 
tionship was obtained between urban NO* emissions derived from 
OMI columns with EDGAR bottom-up inventory estimates. In the 
Saudi Arabian city of Riyadh a strong dependency of the N0 2 column 
density on wind direction was observed, with the maximum con- 
centration shifted by ~10— 40 km downwind from the city center. 
Across all the cities, the NO* lifetimes were in the range of 2.3— 6.4 h, 
reflecting the oxidizing capacity of the freshly emitted plumes. 

As indicated earlier, it is possible that U.S. air quality studies 
could be improved by the use of satellite data to augment uncertain 
emission inventories for neighboring countries, of which Mexico is 
the prime example. Some satellite studies have already shed light 
on emission estimates in Mexico. The emissions inventory of NO* 
was evaluated during the Megacity Initiative: Local And Global 
Research Observations (M1LAGRO/INTEX-B) campaign with a 
combination of in situ aircraft measurements, satellite retrievals, 
and numerical modeling (Boersma et al., 2008a, 2008b; Molina 
et al., 2010). OMI top-down emission estimates have been 
compared with the official Mexican National Emissions Inventory 
of 1999. Large sources of S0 2 were identified on either side of 
Mexico City during the same campaign (de Foy et al., 2009), using 
the oversampling technique described in Section 4.2. Source 
strengths were compared with surface measurements using 
mesoscale simulations, improving the confidence in prior emis- 
sions estimates. Shim et al. (2009) used 0 3 and CO data from TES to 
study pollution outflow from Mexico City, also as part of the 


MILAGRO/INTEX-B campaign. Retrieval issues prevented measure- 
ments directly over the metropolitan area. Generally, TES 
enhancement ratios of 0 3 to CO were overestimated compared to in 
situ measurements, prompting a call for further examination of the 
data sets. 

4.6. Natural point sources 

Volcanoes are the only natural point sources of large magnitude, 
but they can be so large that they can dwarf the emissions of S0 2 , 
PM, and other gases from anthropogenic point sources during 
eruptive and even degassing periods. Sizeable eruptive events can 
release several Tg of S0 2 in a day or two, as compared with about 
1 Tg per year from large anthropogenic point sources like the 
copper and nickel smelters at Norilsk (see, e.g., Walter et al., 2012). 
The history of observations of volcanic emissions from space has 
been reviewed by Thomas and Watson (2010), and the reader is 
referred to that source for more information. One clear advantage 
that satellites have over ground measurements is the ability to 
monitor emissions from a source that may be very remote and 
inaccessible, as well as being expensive and hazardous for humans 
to monitor in situ. Because the volcanic source strength is great, the 
S0 2 signal is readily detectable by several instruments and more 
easily converted to an estimate of emissions than the S0 2 released 
by anthropogenic sources. 

Satellite observations of S0 2 releases from volcanoes trace back 
to the late 1970s and are perhaps best exemplified by the obser- 
vation by TOMS of the eruption of the El Chichon volcano in 
southern Mexico in 1982 (Krueger, 1983). It was estimated from the 
UV absorption properties of the released cloud that ~3.3 Tg of S0 2 
were injected into the stratosphere over a two-day period. One of 
the first comprehensive reports of satellite observations of S0 2 
emissions from volcanoes was by Khokhar et al. (2005), who 
studied both eruptions and degassing at about 20 volcanoes around 
the world during the period 1996—2002, using GOME data. Despite 
the coarse spatial resolution of GOME and the relatively large un- 
certainties for individual observations, it was felt at that time that 
there was potential for being able to derive rough estimates of S0 2 
emissions from satellite retrievals. Uncertainties could be reduced 
if additional information on plume height, cloud cover, aerosol 
loading, and S0 2 lifetime were available. 

As the next generation of satellites came into operation, having 
greater spatial resolution, it became possible to develop additional 
information about volcanic releases of S0 2 . For example, Loyola 
et al. (2008) studied four volcanic eruptions in Central and South 
America with a combination of GOME, SCIAMACHY, and GOME-2 
instruments. By combining the retrievals with a trajectory model, 
they were able to determine the effective plume heights. They also 
proposed the development of a near-real-time S0 2 data delivery 
service to warn of aviation hazards arising from volcanic ash clouds. 
Using OMI and other retrievals, Carn and Prata (2010) reported on 
the numerous episodes of explosive S0 2 emissions on the Carib- 
bean island of Montserrat during the period 1995—2009 and 
particularly in 2003 and 2006. The largest single emission was 
~2 Tg in 2006 that was injected into the lower stratosphere and 
tracked halfway around the world. They suggested that space- 
based routine monitoring of gaseous releases from volcanoes 
might be capable of warning of impending eruptions. 

The August 2008 eruption of the Kasatochi volcano on the 
Aleutian Islands of Alaska, released similar amounts of S0 2 , esti- 
mated at ~ 1.5— 2.5 Tg. Krotkov et al. (2010) combined OMI satellite 
data with a trajectory model to analyze the dispersion of the S0 2 
cloud. The OMI data were able to constrain the S0 2 plume height 
and eruption times and initialize the transport model. The OMI S0 2 
burden was integrated over the Northern Hemisphere and 
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extrapolated back to yield an initial erupted mass of ~2.2 Tg, twice 
as large as the measured mass on the day. In addition, instruments 
on the CAR1BIC flying observatory detected remnants of the Kasa- 
tochi SO2 plume in central Europe, and comparisons with GOME-2 
satellite observations were possible within just a few hours, 
yielding a robust correlation (Heue et al., 2010). Clarisse et al. 
(2011a) were able to observe H2S in the Kasatochi plume from 
IASI observations, though H2S is spectroscopically difficult to detect 
and it oxidizes to S0 2 in the atmosphere within a few days. 

The eruption of the Icelandic volcano Eyjafjallajbkull in 2010 
enabled Carboni et al. (2012) to develop a new scheme to utilize 
IASI retrievals to track the SO2 plume in the presence of cloud and 
volcanic ash. Rix et al. (2012) were able to study S0 2 and BrO col- 
umns in the eruption from this volcano using GOME-2 data, as well 
as to determine the S0 2 plume height. Finally, Martinez-Alonso 
et al. (2012) reported the first observation of CO emissions from a 
volcano using MOPITT data. 

4.7. Natural area sources 

There are several types of natural area sources that are 
important to characterize from the perspectives of assessing the 
contribution of anthropogenic emissions to total emissions, un- 
derstanding atmospheric chemistry and the production of sec- 
ondary species, and assembling complete emissions datasets to run 
CTMs and other atmospheric models. Such source types include 
open biomass burning, soils, lightning, and biogenic sources. 
Biogenic NMVOC were described in Section 3.5 and are not dis- 
cussed further here. 

Satellite estimates of open biomass burning have become 
routine, and their resulting emissions are now reported on a regular 
basis, e.g., the Global Fire Emissions Database, GFED (van der Werf 
et al., 2010), the Fire INventory from NCAR, FINN (Wiedinmyer 
et al., 2011), and the Wildland Fire Emission Inventory, WFEI 
(Urbanski et al., 2011). The literature related to biomass burning 
detection from space is vast and beyond the scope of this article. 
Given the large uncertainty in bottom-up estimates of CO emis- 
sions, inversion analyses using CO satellite observations offer 
considerable promise for constraining this important source (e.g., 
Petron et al., 2004) and separating anthropogenic and biomass 
burning components. 

Lightning NO* has been studied by a number of research groups 
using satellite products. A typical approach to this problem is to 
measure the SCD of lightning-related NO* and then develop a 
formalism to relate SCD to VCD. Cloud-resolving models are useful 
to calculate the appropriate AMF. A further complication is that 
most of the gas exists in the upper troposphere in the form of NO, 
which is not directly measurable. Beirie et al. (2009) identified the 
uncertainties and sensitivities involved in using satellite observa- 
tions to study lightning NO*. The analytical process involves 
generating an estimate of production efficiency (number of mole- 
cules produced per flash) and then multiplying that by an estimate 
of the number of flashes in a given time period over a given area to 
yield an estimate of the NO* produced. 

Perhaps the most widely accepted estimate of annual NO* pro- 
duced from lightning worldwide developed by a method that does 
not involve satellite observations is ~5 Tg N in the range of 2— 
8 Tg N (Schumann and Huntrieser, 2007). There have been a 
number of satellite-based studies that have attempted to constrain 
the global magnitude of lightning NO* generated, and we touch on 
just a few of them here. Though uncertainties are still large, central 
estimates agree quite well. More than a decade ago, Bond et al. 
(2002) used satellite imagery from the Lightning Image Sensor 
(LIS) on board the Tropical Rainfall Measuring Mission (TRMM) 
satellite to obtain an emission estimate of 6.3 Tg N yr _1 in the 


tropics. Beirie et al. (2004c) combined the LIS data on lightning 
activity with GOME NO2 measurements for a statistical case study 
over central Australia, where other NO* sources are very few. 
Extrapolation of the results yielded a global production estimate of 
2.8 (0.8—14) Tg N yr _1 Boersma et al. (2005) also used GOME 
measurements to estimate an annual global NO* production rate of 
1.1 -6.4 TgN in 1997. 

Martin et al. (2007) used observations from four satellite in- 
struments to constrain the estimate of NO production by lightning: 
NO2 columns from SC1AMACHY, O3 columns from OMI and MLS, 
and upper tropospheric HNO3 from ACE-FTS. Supplementing these 
observations with GEOS-Chem model predictions of the times and 
places of lightning occurrences, a global estimate of 6 ± 2 Tg N yr _1 
was obtained. When Beirie et al. (2010) studied SCIAMACHY mea- 
surements of lightning NO* over the eastern U.S., they often found 
poor spatial correlation between fresh flashes and enhanced col- 
umn densities, inferring that aged lightning NO* or anthropogenic 
NO* was interfering. Production efficiency was often lower than 
expected, equivalent to a global source strength of ~1 Tg N yr 1 
Bucsela et al. (2010) were able to develop customized OMI re- 
trievals for lightning NO* during the NASA TC 4 campaign, which 
were supported by aircraft observations. Combining CMAQ with 
OMI observations, an estimate of 500 mol NO per flash provided 
reasonable (-5 to +13%) model/measurement agreement for the 
mean summertime tropospheric N0 2 column (Allen et al., 2012). 

Soil NO* emissions have also been estimated from satellite ob- 
servations. It is difficult to estimate soil NO* emissions by bottom- 
up methods, due to high spatial and temporal variability, as well as 
the sensitivity of the emission rates to climatic variables, soil 
characteristics, and anthropogenic influences such as cultivation 
and fertilizer application. Soil NO* emissions might therefore be a 
prime candidate for measurement from space, where spatial 
coverage is high and daily observation is possible. As with lightning, 
NO is the primary emitted species and the rate of oxidation to other 
nitrogen species is an issue. Jaegle et al. (2004) first demonstrated 
the ability of satellites to observe soil NO* emission processes. 
Bertram et al. (2005) studied an agricultural area of two million 
hectares in Montana using SCIAMACHY data and observed pulses of 
NO* associated with application of fertilizer and precipitation. They 
were able to use the satellite data to refine mechanistic parameters 
in a conventional model of soil NO*. Hudman et al. (2010) used OMI 
data to study soil NO* emissions in the Great Plains. They also used 
a case study of a region of soybean and corn production in South 
Dakota to demonstrate rain-induced pulses of NO*, as well as as- 
sociations with fertilizer application. Large interannual variation in 
the emission rate was observed during the period of measurement 
(2005—2008), and differences between model simulations and 
satellite retrievals were noted. Overall, the soil NO* emissions 
resulted in an increase of ~3 ppbv in 8-h O3 concentrations over 
the agricultural Great Plains region. Subsequently, Hudman et al. 
(2012) were able to develop a mechanistic model of global 
soil NO emissions using OMI satellite data to constrain the 
parameterization. 

4.8. Field campaigns: observing transport and transformation 

Satellite observations are becoming an increasingly integral part 
of field campaigns aimed at understanding air quality and source 
contributions. Although evaluating and improving space-based 
measurement of tropospheric composition is a nearly universal 
goal of such field experiments, satellite products can be part of the 
planning process and can contribute to analysis of the results. 
Consider 1NTEX-A (Intercontinental Chemical Transport Experi- 
ment-North America) with the objective of studying pollution 
plumes over the U.S. and export from North America, and INTEX-B 
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(Intercontinental Chemical Transport Experiment— Phase B) with 
the objective of quantifying the outflow of pollution from Mexico 
and East Asia. These projects integrated satellite products along 
with model forecasts into the flight planning process. In TexAQS II 
(The Second Texas Air Quality Study) satellites provided the 
regional context for urban measurements of trace gases (McMillan 
et al., 2005) and aerosols (Winker et al., 2007). 

In several recent, large-scale field campaigns, satellite observa- 
tions have been a core component. EAST-AIRE (East Asian Study of 
Tropospheric Aerosols: An International Regional Experiment), 
RAMMPP (Regional Atmospheric Measurement, Modeling and 
Prediction Program) and DISCOVER-AQ (Deriving Information on 
Surface conditions from Column and Vertically Resolved Observa- 
tions Relevant to Air Quality) offer illustrative examples. For these 
projects, satellite observations were incorporated into the experi- 
mental plan as one of the tools to understand atmospheric pro- 
cesses. EAST-AIRE (and its successor EAST-AIRC) were focused on 
aerosol emissions, transport, distributions, optical and chemical 
properties, and their subsequent impact on direct radiative forcing, 
air quality, ecosystems, and regional climate, as summarized in Li 
et al. (2007, 2011). A 2008 experiment employed a combination 
of aircraft and satellite observations to constrain the S budget over 
China and to estimate the lifetime and export (He et al., 2012). The 
combination of aircraft and satellite observations allowed the 
plume to be tracked for several days, resulting in an estimate of the 
gas-to-particle conversion rate (Dickerson et al., 2007; Krotkov 
et al., 2008; Li et al., 2010a). 


RAMMPP is a long-standing program to provide policy-relevant 
science to the nonattainment areas of the Mid-Atlantic states 
(Castellanos et al., 2011: Hains et al„ 2008; Lee et al., 2011a; 
Loughner et al., 2011; Marufu et al., 2004; Taubman et al., 2006). 
The rate of emission of S0 2 over the U.S. is reasonably well known, 
because many large sources are monitored with CEMS. The lifetime 
of the S0 2 , however, depends on concentrations of the oxidants OH 
and H 2 0 2 as well as cloud cover, variables that are hard to model 
accurately with CMAQ (Loughner et al., 2011). Because S0 2 is the 
major precursor of sulfate particles, the rate of conversion must be 
understood in compliance modeling for PM 2 5 . Both in situ and 
remotely sensed data have been used to estimate the overall rate of 
loss of S0 2 from the atmosphere and by implication the rate of 
formation of sulfate (Lee et al., 2011a). Both RAMMPP and INTEX-A 
(Hennigan et al„ 2006) show a nearly exponential decrease in 
mixing ratio with altitude with the bulk of the S0 2 located below 
2000 m. A mass balance approach indicates an S0 2 lifetime of 17 
(±7) hr for summer conditions. Lee et al. (2011a) expanded these 
observations to the global scale using GEOS-Chem and column 
contents for OMI and SCIAMACHY. The results (Fig. 8) show the 
seasonal cycle and latitudinal dependence expected for photo- 
chemical oxidation processes, with maximal lifetimes in the winter 
and at high latitudes. 

A major challenge for Earth-observing satellites is to distinguish 
between pollution high in the atmosphere and pollution near the 
surface that has a direct connection to emissions and a direct 
relevance to human health (Lin and McElroy, 2010). In summer 




Fig. 8. Modeled and measured boundary-layer S0 2 tropospheric lifetimes over the eastern U.S. and global lifetimes of S0 2 using a combination of numerical modeling as well as in 
situ and remote observations (adapted from Lee et at, 2011a). This material is reproduced with permission of John Wiley & Sons, Inc. 
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2011, the DISCOVER-AQ program began a multi-year airborne field 
campaign to tackle this challenge. The first experiments were 
performed over the urban area encompassing Baltimore, MD and 
Washington, DC. Results indicate that there is a strong correlation 
between column content and surface concentration, but 0 3 located 
above the PBL can interfere, especially on days with poor air quality. 
The main precursor of 0 3 formation over the eastern U.S. is N0 2 . 
Preliminary results indicate that the bulk of N0 2 lies near the sur- 
face, both in the morning and afternoon. If N0 2 is generally 
confined to the surface, then satellite observations should be ex- 
pected to correlate well with surface observations. The similarities 
of NO* both upwind of the Baltimore/Washington area (morning 
measurements) and downwind (afternoon measurements) indicate 
that O3 is formed over a broad area. This shows that emissions 
controls can yield benefits hundreds of kilometers from the source. 
Results from OMI (Fig. 9) indicate high concentrations of N0 2 near 
sources such as New York City, but substantial concentrations over 
rural areas as well. CMAQ captures the major source areas seen by 
OMI. The ratio of urban to rural N0 2 column content, however, is 
higher in CMAQ than in OMI or the in situ data (Allen et al„ 2012). 
Current work involves verifying these differences and exploring 
how the dynamics and chemistry of the CMAQ model could be 
improved. Satellites can provide coverage unavailable from even 
the most ambitious in situ measurements and when compared 
appropriately to numerical simulations can provide guidance on 
model improvement. Because plans for compliance with air quality 
standards rely on simulated emission reductions, such models must 
be able to reliably reproduce the temporal and spatial character- 
istics of pollution events. 

5. The National Emissions Inventory (NEI) 

It is clear from the content of Section 4 that there has been a 
tremendous amount of work undertaken in the past decade to 
interpret satellite retrievals with the purpose of gaining new in- 
formation about the spatial and temporal aspects of emissions of a 
variety of species from many types of sources. The value that these 
studies have depends partly on what the perceived needs are for 
the regulatory community and the degree of confidence that we 
have in bottom-up inventory methods. At the end of the day, the 
capability of the satellite instruments must be weighed against the 
established inventory methods currently in use. Only then can we 
judge the merits of focusing on particular sources and species and 


deciding what resources should be invested in pushing the appli- 
cation of satellite methods into the emissions estimating arena. 

The U.S. National Emissions Inventory (NEI) is a national 
compilation of emissions sources collected from state, local, and 
tribal (S/L/T) air agencies specifically for the NEI, as well as from 
other EPA emissions programs, including the Toxics Release In- 
ventory (TRI), emissions trading programs such as the Acid Rain 
Program, and data collected as part of EPA regulatory development 
for reducing emissions of air toxics. The NEI development efforts 
include assembling component datasets, blending data from mul- 
tiple sources, and performing quality assurance steps that further 
enhance and augment the compiled data. The NEI is produced by 
EPA on a triennial cycle, with the most recently released NEI for the 
calendar year 2008. 

The species included in the NEI are those related to imple- 
mentation of the National Ambient Air Quality Standards 
(NAAQS) — known as criteria air pollutants (CAPs), other precursors 
of ozone and PM, as well as hazardous air pollutants (HAPs) asso- 
ciated with EPA's Air Toxics Program. The CAPs and their precursors 
include lead (Pb), CO, NO*, VOC, S0 2 , NH3, PM of diameter 2.5 pm or 
less (PM 2 5 ), and PM of diameter 10 pm or less (PMi 0 ). The HAPs 
include the 187 remaining HAPs from the original 188 listed in 
Section 112(b) of the 1990 Clean Air Act Amendments. Key HAP 
emissions species include mercury (Hg), hydrochloric acid (HC1) 
and other acid gases, heavy metals such as nickel and cadmium, 
and hazardous organic compounds such as benzene, formaldehyde, 
and acetaldehyde. 

While the NAAQS program is the basis on which EPA collects 
CAP emissions from the state, local, and tribal air agencies, it does 
not require collection of HAP emissions. For this reason, the HAP 
reporting requirements are voluntary. Nevertheless, the HAP 
emissions are an essential part of the NEI program. These emissions 
estimates allow EPA to assess progress in meeting HAP reduction 
goals described in the Clean Air Act amendments of 1990. These 
reductions seek to reduce the negative impacts to people of HAP 
emissions in the environment, and the NEI allows EPA to assess 
how much emissions have been reduced since 1990. With the 
exception of a few organic compounds like formaldehyde, HAP 
species are not at present able to be detected by satellite in- 
struments, but that may change in the future. 

The Air Emissions Reporting Rule (AERR) is the rule that requires 
states to submit estimates of CAP emissions and provides the 
framework for voluntary submission of HAP emissions. The AERR 
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Fig. 9. N0 2 columns for the month of July 2011 over the Mid-Atlantic states during DISCOVER-AQ. The CMAQ model matches urban observations well but underestimates rural N0 2 . 
Increasing the photolysis rate of alkyl nitrates in CB05 improves model-measurement agreement. 
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has emissions thresholds that determine whether an S/L/T agency 
must report stationary source emissions as point sources or 
whether the emissions can be lumped together into a county total 
as nonpoint sources. Based on the AERR requirements, S/L/T 
agencies submit emissions of point, nonpoint, on-road mobile, 
nonroad mobile, and fires emissions sources. These submissions are 
sent to the Emission Inventory System (EIS) that EPA then uses to 
review and assemble the data from the S/L/T agencies. 

While the basis of the NEI is an S/L/T agency reporting program, 
EPA compiles data for use by these agencies, including data that are 
derived, or could be derived, from satellite-based information. 
While initial resource investments are generally needed for 
updating emissions methods to use satellite data, there can be long- 
term benefits of the resultant improved data, such as making better 
regulatory or other decisions. 

Some uses of satellite data already exist in the NEI. A satellite- 
based emissions method is the use of fire detect data from satel- 
lites in the SMARTFIRE version 2 system as described in EPA 
(2013a). While the satellite data are useful for compiling reason- 
able location and temporal information, the data are best used for 
emissions estimation when accompanied by ground-based obser- 
vational information to augment the spatial extent of the fires as 
well as the fuel consumption estimates derived from other models. 
In addition, satellite images are routinely used by S/L/T agencies 
and EPA to verify location information for point-source facilities. 

Other satellite data may be able to be used in the future for the 
NEI to compile other types of activity data. For example, satellite 
images could conceivably be used to develop or augment traffic 
count data, the size and location of animal feedlots and waste la- 
goons, the locations of residential and non-residential construction, 
and the characterization of “urban” or “rural” that is a part of some 
methods such as residential wood combustion. Depending on the 
frequency of availability, satellite images could also be used to help 
assess the timing of activities such as agricultural tilling and asphalt 
paving, where the timing relative to meteorological events can 
determine whether the activity actually occurs and what the 
emissions will be. 

A summary of the construction and the components of the 2008 
NEI by sector, key species, and sources of data is presented in 
Table 2, with notes regarding uncertainties associated with the 
various sectors. The strengths and weaknesses of the NEI vary from 
year to year as changes are made to the creation processes and as 
resources for the program vary. It is impossible to generalize about 
the strengths and weaknesses since data are compiled from so 
many sources. Since the most recent NEI was released for calendar 
year 2008, this summary reflects only on that inventory. EPA esti- 
mates emissions from many source categories and then uses those 
data when states do not submit it; thus, the quality and origin of the 
data can vary from state to state. As a general matter, the primary 
sources of uncertainty are: ( 1 ) uncertainty in emission rates; ( 2 ) 
uncertainty in the activities generating emissions; and (3) the 
temporal and spatial variability of both emission rates and 
activities. 

6. Promising applications and future needs 

Organizations that develop emission inventories can profit from 
the vast amount of spatial and temporal data delivered from sat- 
ellites; however, much of this opportunity has yet to be realized. 
One challenge is that emission inventories consist of a series of 
data-intensive calculations, while the results from most inverse 
modeling studies consist of two-dimensional maps of emission 
adjustments. Progress on the use of satellite data to inform emis- 
sion inventories depends on connecting these emission adjust- 
ments to the real building blocks of emission inventories, namely 


activity levels and emission factors. This is possible with further 
collaboration between emission inventory developers and the 
satellite research communities. In this final section we review some 
ways in which the gap between promise and fulfillment of satellite 
applications for emissions estimation might be narrowed. 

Identification of the practical problems faced routinely by 
emission inventory developers can be used to guide further inter- 
action between the communities. For example, the acquisition of 
temporal information from satellites to help adjust day-of-week 
emission allocations that are needed for Eulerian grid models 
would be beneficial for many types of area sources, including 
onroad mobile, nonroad mobile, and shipping. A number of un- 
conventional area sources, e.g., evaporation ponds from oil and gas 
development activities, are desperately short of data. Often there is 
industry and state resistance to monitoring such sources. If the 
detection limit for such emissions could be improved, there is great 
potential to enhance the conventional bottom-up inventory 
methods. State regulators have expressed a need for help in iden- 
tifying small but significant hotspots of emissions of species such as 
SO 2 . Monitors are sparse in many areas and expensive to maintain. 
This means that some sources may elude detection by the ground 
monitoring network. To be able to get uniform coverage over a very 
wide area could be extremely valuable and cost-effective. 

Ammonia is difficult to measure in situ and measurements are 
sparse; but reduced N species are becoming an increasingly 
important component of reactive N deposition and are steadily 
leading to adverse environmental effects. Detection of agricultural 
emissions of NH 3 from space would be an ideal opportunity to 
improve inventories and identify hotspots. Observations could be 
used to help ground-truth the results that will come from the 
updated NH 3 emission rate information provided by the National 
Air Emissions Monitoring Study (NAEMS). And there is a continued 
need in the regulatory community for satellite products that allow 
estimation of biomass burning activities using software such as 
SMARTFIRE to support a nationally consistent approach that brings 
together the best of the satellite and on-the-ground information for 
a robust inventory. 

Using ratios of trace gases as an indication of emission amounts 
and characteristics also shows great promise, but these require co- 
location of observations. For example, biomass burning and other 
inherently inefficient combustion processes have elevated CO/C0 2 
and CH 4 /CO 2 ratios, while modern, efficient power plants have very 
low ratios. Unregulated or under-regulated internal combustion 
engines can produce elevated NO 2 /CO 2 ratios. Such information can 
help to improve emission inventories and help policy makers target 
the most important sources. 

The temporal and spatial characteristics of natural emitting 
sources are exceedingly difficult to obtain by bottom-up methods 
and are usually highly uncertain; however, they are important to 
know because they can represent a sizeable fraction of total pri- 
mary emissions and are known to play a highly significant role in 
atmospheric chemistry. They also contribute in important ways to 
the ambient concentrations of many primary and secondary species 
in the atmosphere. Satellite measurements have so far proven 
valuable in application to a number of types of natural sources: 
biogenic NMVOC, volcanic S0 2 and PM, wetlands CH 4 , biomass 
burning CO and PM, soil and lightning NO x , etc. To the extent that 
information about these emitting sources is necessary for modeling 
their contributions to ambient concentrations, they are a necessary 
part of the information portfolios of air quality managers, and the 
benefit of space-based observations has not yet been fully realized. 

Satellite detection of aerosol speciation is an important future 
need from the public health perspective. A clear understanding of 
the spatial trend and temporal variation of PM species is crucial for 
U.S. EPA to design cost-effective emission control policies, as well as 
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Table 2 

Summary of the 2008 NEI by sector, key species, and sources of data. 


Sector 


Species Sources of data 


Agriculture 


NH 3 , PMio, PM 2.5 


Dust — Construction, paved, PMio, PM2.5 

and unpaved road dust 

Fires CO, VOC, PM 2 . 5 , 

PMio, HAPs 


Fuel Combustion — Electric 
Generation 


NO x , S0 2 , PM 2 . 5i 
H g, HAPs 


Fuel Combustion — Industrial, NO*. S0 2 , PM 2 5 , 
Commercial, and Hg, other HAPs 

Institutional 


Fuel Combustion — residential CO, NO x , PM 2 5l 
PMio, HAPs 


Industrial Processes 


NO x , VOC, S0 2 , 
PM 2 .s, HAPs 


Onroad mobile 


CO, NO x , VOC, 
HAPs 


Nonroad mobile equipment CO, NO x , VOC, 
HAPs 

Nonroad mobile aircraft, Pb, (aircraft), 

railroad, and marine CO, NO x , S0 2 


Solvents 


VOC, HAPs 


Nonpoint with minor point-source contributions by a few states who submit farms as point sources. Blend of 
state and EPA data, with about 75% EPA estimates. NH 3 emission rate improvements pending ongoing NAEMS 
study. Ambient data suggest dust emissions are too high. In addition to emission rates, sources of uncertainty 
include: (1) the timing of animal waste operations, which impact the emission rates in combination with 
meteorological influences, (2) the location of animal operations, (3) the timing of tilling operations, (4) which 
fields are tilled in a given year, (5) the timing of fertilizer application, and (6) which fields are fertilized in a 
given year. 

Nonpoint with very minor point-source contributions from industrial sites. Mostly EPA estimates, generally 
about 70%. Meteorological adjustments applied at monthly resolution to consider moisture. Ambient data 
suggest dust emissions are too high. Paved roads emission rates updated in 2008 NEI based on latest data. 
Satellite-based fires estimates from SMARTFIRE version 2, including allocation of all fires to wildfires, 
prescribed burning, and agricultural fires. Emissions are both spatially and temporally resolved by day. 
Emissions rates may be out of date, but documented. HAP emission rates do not vary spatially. Satellite-based 
methods are known to both overestimate and underestimate areas burned. Sources of uncertainty include the 
assessment of available fuel, fuel consumption, quantification of smoldering, extent of a burn, and the type 
of the fire (a wildfire or a prescribed fire). 

Data are largely based on data from continuous emissions monitors (especially in the Eastern US), which 
measure S0 2 and NO x (and in few cases Hg) emissions. These data and heat input data are readily available 
from the EPA and are used by states whose data submissions make up nearly 100% of the inventory for 
NO x , S0 2 . Other species from both EPA and states, with EPA data from the Mercury and Air Toxics Standard 
being used for most EGUs for Hg, other metals, and acid gases. Sources of uncertainty include emissions from 
sources that do not have continuous emissions monitors. 

Includes a broad mixture of point and nonpoint data, state and EPA data, which makes generalizations about 
the data sources difficult. While nonpoint estimates generated by EPA are based on fuel consumption data 
from the Department of Energy, fuel consumed by sources reported by states as point sources needs to be 
subtracted from the nonpoint calculations. While states were responsible for this step in 2008 NEI, we do 
not know if this was always done and suspect it was not in some cases. Some states may have double 
counted, while other states may have left out nonpoint emissions. Emission rates are outdated in some cases 
(particularly for PM 2 .s) with new methods now available (but not for 2008 NEI) for condensible PM 2 .5 
emissions. Some ( ~ 1 5%) HAP emissions from industry reporting to the EPA Toxics Release Inventory. 

Data are largely a nonpoint estimate that can be reasonably estimated based on fuel use and emission 
factors. Natural gas heating important for NO x and is about 80% reported by states and 20% EPA estimates. 
Wood heating is important for PM 2 . 5 and is about 55% reported by states. Urban area wood combustion 
estimates have been called into question with overestimates apparent as compared to ambient data. 
Methodologies and default emission factors revised by state/EPA and improved for 2008. 

Almost all data provided by states as point sources. Quality varies greatly depending on the source of the 
information provided to EPA, which is generally not possible to assess because the origin of the data is usually 
not provided. For HAPs, EPA augments with TRI data submitted by industry to EPA, but this is 10—30% 
depending on the industry. 

Estimated by the onroad emissions model called the Motor Vehicle Emissions Simulator (MOVES), based on 
modeled hourly temperatures using a 12-km gridded modeling resolution. California emissions provided by 
that state. For MOVES-based emissions, data are computed hourly and summed up for the year. Many states 
did not have state-specific inputs available because of the newness of the model. Activity data are based on 
Federal Highway traffic count data at the state level and apportioned to counties, road types, and vehicles. 
Estimated by the NONROAD model, except in California. Some known problems with double counting of 
emissions in a few states, though this is atypical for this sector. EPA keeps emission rates updated with the 
latest emission rates and equipment populations and activity estimates. 

Primarily EPA data for key species, with a relatively high contribution of state-submitted data for NO x from 
aircraft of 15%. EPA data are estimated at a national resolution and apportioned to counties, rail lines, ports, 
shipping lanes, and airports using methods developed in collaboration with many states. Emission factors 
are updated regularly by EPA, but apportioning approaches for category 1&2 marine vessels, which have in 
some cases overestimated emissions in ports in some demonstrated cases as compared to ambient data. 

100% state data for most solvent categories, with EPA not augmenting these sectors due to overlap between 
point and nonpoint emissions that the states are expected to resolve. Have confirmed cases where states did 
not submit VOC HAPs but did submit VOC and EPA did not compute or add HAP VOC. Methodologies and 
default emission factors revised by state/EPA and improved for 2008. Emission factors for HAPs out of date 
and not consistent with new product formulations and technology changes. Methods updated for 2008 with 
activity data on solvent use used from Freedonia, though not for calendar year 2008. 


to update NAAQS species and HAPS that have different primary 
emission sources and different toxicity levels. EPA’s PM mass 
monitoring network has ~ 1200 sites, covering less than 25% of the 
3100 counties in the continental U.S. In addition, EPA’s Chemical 
Speciation Network and IMPROVE network together have only 
~500 sites. This number is unlikely to increase given current 
budgetary constraints. Satellite instruments such as MISR and OMI 
have shown promising signs of ability to differentiate boundary- 
layer fine PM components such as sulfate and black carbon. How- 
ever, new algorithms or new sensors are needed to achieve the 
level of retrieval accuracy and spatial resolution that would be 


useful for air quality applications. For example, given the lower 
tropospheric lifetime of fine PM (a couple of days to a week), a 
satellite must be able to achieve a return frequency of 1—2 days in 
order to study air pollution episodes related to wildfires, dust 
storms, and heavy haze. The large spatial gradient of black carbon 
requires the sensor to have a spatial resolution of 1 —2 km. 

As data from the current generation of satellites are used more 
and more in tropospheric pollution applications, it is becoming 
apparent that more intensive validation efforts are required. Two 
major factors limiting our confidence in inversion studies are 
uncharacterized model errors, which are addressed in Section 4.1, 
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and systematic errors in the satellite retrievals. Typically, we only 
assume the presence of a 0th order systematic error (bias) in the 
satellite products. However, satellite retrieval errors can be corre- 
lated with other geophysical quantities, such as temperature, water 
vapor, surface reflectance, and cloud properties. To properly di- 
agnose these potential systematic errors, we need coincident ob- 
servations of both the species retrieved and any geophysical 
quantities that could cause systematic errors in the retrieval. At a 
fundamental level, when we assimilate observations in an inverse 
modeling framework, we want to be confident that the noise in the 
observations is truly random. 

All air quality observations from space so far have been from low 
Earth orbit (LEO). These provide return times for observation of at 
most once per day at a specific time of day. Because of clouds and 
retrieval noise, useful data may therefore be obtainable only with 
weekly or monthly temporal resolution. Observation from geosta- 
tionary Earth orbit (GEO) would be of considerable value for air 
quality investigations by providing continuous mapping on a con- 
tinental scale with hourly resolution and on the order of 
4 km x 4 km spatial resolution. As mentioned in Section 2, GEO- 
CAPE and TEMPO are planned geostationary missions to observe 
tropospheric air pollution over North America. However, with 
current launch dates for both these missions likely to be 2017— 
2020 or later, acceleration of the schedules would greatly benefit 
the U.S. air quality management community in a number of ways. 
Continuous atmospheric composition from GEO measurements 
would provide for the first time the data density required for a 
comprehensive investigation of pollution events, pollution trans- 
port, and atmospheric processes at different altitudes and across 
fine spatial and temporal scales. In particular, such data would 
permit inversion of emissions at the high spatial and temporal 
resolution necessary to characterize diurnal variations. Observa- 
tional errors for the inversion could be effectively decreased 
through the averaging of a large number of retrievals collected for a 
given scene during integration times of an hour or less. The high 
spatial and temporal resolution achievable in inversions would 
lend insights into the causes of bottom-up and top-down in- 
consistencies and potential errors in current emission inventories. 
They would also provide a unique ability to constrain transient and 
episodic sources such as wildfires and industrial chemical releases. 

To quantify and separate the influences on chemical distributions 
of varying emissions, photochemistry, and dynamical processes, 
it is desirable to provide improved vertical profile information 
from satellite observations, especially in the lowermost troposphere. 
This has been demonstrated from LEO for CO (Worden et al., 2010) 
and O3 (Fu et al., 2013) through the use of multispectral retrieval 
techniques. Similar measurements are proposed for the GEO-CAPE 
mission, and an independent estimate of concentrations in the 
lowest 2 km should significantly improve the capability to constrain 
surface emissions. 

In his eulogy for the loss of SCIAMACHY (IGAC News, August 
2012), John Burrows reminded us that “we cannot manage what 
we can’t measure.” Reflecting on this adage is a useful way to 
conclude this article. U.S. air quality managers must continually 
manage air quality, even in the face of measurement shortcomings. 
Emission factors are measured for a few sources and then extrap- 
olated to the array of real-world devices; emission rates are 
measured at the stack for those sources deemed large enough to 
warrant it; activity levels are inferred from surveys of the number 
of sources and their expected patterns of operation; and, in a few 
cases, monitors are set up in the field and source emissions deter- 
mined by modeling. The situation is far from ideal, but managers 
have no choice but to accept the methodological limitations and 
acknowledge the levels of uncertainty that result. Perhaps it would 
be more appropriate to say that “what we can’t measure well, we 


cannot manage effectively.” The lure of satellite observations is that, 
given the right conditions, we might one day be able to measure the 
actual emissions from many sources at the same time, all of the 
time. 
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